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I •  INTRO DUCTION 

A.  BACKGROUND 

For  a  few  years  the  Surface  Effect  Ship  (SES)  has  been 
receiving  increasing  attention  in  the  United  States  Navy 
(USN)  ,  and  detailed  studies  have  been  undertaken  since  then. 
The  two  major  categories  of  SES  are  the  Air  Cushion  Vehicles 
(ACV)  and  the  Captured  Air  Bubble  Ships  (CAB) . 

Whereas  the  weight  of  the  ACV  is  entirely  supported  by 
the  pressure  differential  in  the  plenum  chamber  all 
surrounded  by  flexible  seals  between  the  water  surface  and 
the  hull  structure,  the  CAB  uses  rigid  sidewalls  extending 
into  the  water  and  thus  giving  additional  support  to  that  of 
the  captured  air  bubble.  Some  lift  forces  are  made  up  by  the 
flexible  bow  and  stern  seals  as  well  as  aerodynamic  forces 
at  higher  speeds  are  usually  taken  into  account. 

The  latter  type  of  craft  has  become  of  primary  interest 
to  the  USN  for  ocean  going  applications  and  led  to  the 
construction  of  several  test  crafts,  among  those  the 
Aerojet-General  100-A,  the  Bell  Aerospace  Systems  100-B, 
both  of  about  100  tons  displacement,  and  the  XR-3  with 
approximately  3  tons  displacement.  This  study  will  be 
concerned  with  the  XR-3. 

This  ship  was  built  by  the  Naval  Ship  Research  and 
Development  Center  (NSRDC)  in  1966  and  tested  for  several 
years.  In  March  1970,  the  XR-3  was  delivered  to  the  Naval 
Postgraduate  School  (NPS) ,  Monterey,  California,  for  further 
investigations.  After  a  digital  computer  simulation  program 
for  the  100-B  test  craft  was  developed  by  Oceanics, Inc . ,  Leo 
and  Boncal  [Ref.  I]  converted  this  program  to  represent  the 
XR-3  test  craft  as  there  were  substantial  design  differences 
between  those  two  ships. 
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In  December  1974,  Forbes  [Ref.  2]  undertook  a  validation 
of  the  Loads  and  Motion  computer  program  in  calm  water  and 
introduced  some  modifications  in  various  program 
subroutines.  At  the  same  time  Finley  [Ref.  3]  refined  the 
bow  and  stern  seal  subroutines  and  developed  new  fan  maps  in 
order  to  achieve  more  realistic  air  flow  rates  and  thus  an 
inproved  craft  performance. 

B.  OBJECTIVES 

The  purpose  of  this  thesis  is  the  study  of  the  pitch  and 
roll  response  of  the  XK-3  Loads  and  Motions  computer  program 
in  a  comparison  between  the  two  modified  versions  of  Refs.  2 
and  3.  Changes  are  introduced,  where  neccessary,  with  the 
idea  to  improve  the  computed  pitch  and  roll  behaviour.  The 
final  goal  is  to  select  an  optimal  model  from  the  results 
obtained  for  future  studies. 

In  order  to  ease  the  further  discussion  the  computer 
models  involved  are  named  as  follows: 

Program  1  =  computer  simulation  program  used  in  Ref.  1 
Program  2  =  computer  simulation  program  used  in  Ref.  2 
Program  3  =  computer  simulation  program  used  in  Ref.  3 

A  users  manual  for  the  final  program  chosen  is  attached 
in  Appendix  A  and  should  be  read  in  conjunction  with  this 
text  in  order  to  better  understand  the  subroutine  programs 
and  variables  under  discussion. 
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II.  GENERAL_DISCaSSION 

A.  ANALYSIS  CF  MOMENTS 

The   development   of   the  equations  of  motion  considering 

the  craft  to   have   six   degrees   of   freedom   will   not   be 

discussed   here   as   it   is  well  presented  in  Refs.  1  and  4. 

However,   the   forces   and   moments    involved   need    some 

consideration.   In   general   ,   various   elements  contribate 

separate  forces  to  the  different  degrees  of  freedom  and   can 

be  represented  by: 

F  =  F  +  F  +F  +F 

bowl  seal     stern  seal     sidewalls     rudder 

+  F  +  F       +  F  +F 

propulsion     waves     aerodynamic    airbubble 

+  F 

gravitational 

With  respect  to   the   local   reference   system   the   x-   and 

z-components   of  the  individual  forces  determine  the  moments 

about  the  lateral  y-axis  and,   therefore,   the   final   pitch 

angle   at   steady  state.  These  moments  are  calculated  in  the 

various  subroutine  programs  (see  Appendix  A)  and  then   added 

up   to   give  the  net  moment  which  eventually  reduced  to  zero 

in  calm  water  steady-state. 

In  figure  1  these  moments  for  the  XR-3  are  shown  in  a 
schematic  manner  to  indicate  their  predominant  effect  on  the 
pitch  behaviour.  Counterclockwise  acting  moments  are 
considered  positive  by  convention. 

In  a  similar  fashion,  the  moments  around  the  x-axis 
determine  the  roll  angle,  which  -  provided  the  craft  is 
symmetric  about  its  longitudinal  axis  including  weight 
distributions  -  tends  to  be  zero  in  calm  water.  In  order  to 
give  an  impression  of  the  magnitude  ratios  the  rounded 
moment  values  for  a  typical  straight  run  at  20  knots  are 
given  in  table  1.  These  figures  were  obtained  from  program  2 
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FIGURE     1 
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FMSS       =    MOMENT  DUE    TO    STERN    SEAL    FORCES 

FMSW       =    MOMENT  DUE    TO    SIDFWALL    FORCES 

FMRUD    =    MOMENT  DUE    TO    RUDDER    FORCES 

FMP         =    MOMENT  DUE    TO    THRUST    OF    PROPELLER 

FMAED    =    MOMENT  DUE    TO    AERODYNAMIC    FORCES 

FMBUB    =    MOMENT  DUE    TO    LIFT    FORCES    OF     AIP    BUBBLE 

FWAVZ    =    MOMENT  DUE    TO    BUBBLE    WAVE    MAKING    DRAG    FORCES 

FMWAV    =    MOMENT  DUE    TO    WAVE    FORCES 

(     EQUALS    ZERO    IN    CALM    WATER     ) 
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in  its  original  version. (See  also  D.1.a.) 


Table  I 

Summary 

of  Moments 

l£i 

tchl 

(Program 

2,  origina 

1  version) 

Component 

Positive 

Negative 

Name 

Moment 

Moment 

(lb-ft) 

(lb-ft) 

FMBS 

705 



FMSS 

--- 

112 

FMSW 



3827 

FMEOD 



85 

FMP 

1147 



FMASD 

645 

FMBUB 

1805 



FWAVZ 



278 

FMWAV 

0 



Summary 

4302 

4302 

B.  INITIAL  CONDITIONS 

Reference  2  points  out  that  the  choice  of  the  program 
initial  conditions  is  somewhat  critical.  If  they  are  not 
close  to  the  steady  state  values  the  imbalance  in  forces 
generated  within  the  program  will  cause  an  execution  stop 
due  to  an  improper  integration  step  size.  Thus  a  tedious 
method  with  elaborate  curves  families  was  used  in  Ref.  2  to 
develop  the  correct  data  including  a  proper  integration  step 
size. 

References  2  and  3,  unfortunately,  give  no  information 
about  the  effect  of  their  program  changes  on  the  initial 
conditions  nor  are  the  values  used  for  the  simulation  runs 
completely  listed.   Refference   2,   however,   presents   some 
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steady  state  values  for  straight  runs  with  pitch  angles 
considerably  larger  than  those  reported  by  01mstead[  F.ef .  5] 
which  come  closer  to  the  values  observed  with  the  present 
craft  configuration. 

As  further  data  were  not  readily  available,  the  initial 
conditions  given  in  Ref.  1  were  chosen  and  are  repeated  in 
table  II.  The  first  trial  runs  with  programs  2  and  3  shoved 
in  both  cases  that  the  disturbances  introduced  were  well 
within  the  limits  of  the  programs'  capability.  Steady  state 
was  reached  after  an  avarage  of  less  than  ten  seconds 
simulated  run  time  of  the  craft. 

Table  II 

Initial  Conditions 
Speed     Pitch    Draft   Plenum     Thrust   Rudder 

9J!3;L§         ££§ssure 
(knots)   (deg)    (in)      (psf)     (lb)    (deg) 

10.0  1.72  6.60  24.86  504.0  0.0 

12.5  1.31  6.50  24.86  420.0  0.0 

15.0  0.94  6.35  24.86  392.0  0.0 

17.5  0.60  6.15  24.86  398.0  0.0 

20.0  0.29  5.98  24.86  424.0  0.0 

22.5  0.13  5.78  24.86  464.0  0.0 

25.0  0.09  5.56  24.86  512.0  0.0 

27.5  0.06  5.30  24.86  568.0  0.0 

30.0  0.05  3.02  24.86  626.0  0.0 

C.     CRAFT    EATA 

The  following  data  were  used  throughout  the  entire  study 
unless  otherwise  specified.  The  numbers  in  parentheses 
indicate  the  block  numbers  under  which  the  values  are  read 
into  the  program  (subroutine  INCON) .  These  names  are  used  in 
this  study  without  further  explanation  which  can   be   looked 
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up  in  Appendix  A. 


XS  (2) 

— 

10.05 

feet 

forward  of 

transom 

XLT0T(3) 

= 

24.7 

feet 

XLSW  (4) 

= 

21.9 

feet 

XSSI  (5) 

= 

3.  97 

feet 

forward  of 

transom 

XLF  (5) 

= 

4.06 

feet 

XBSI  (6) 

= 

2  3.44 

feet 

forward  of 

transom 

XBF  (6) 

= 

4.  06 

feet 

XL3W(7) 

= 

20.0 

feet 

XPKV  (7) 

= 

17.2 

feet 

forward  of 

transom 

XL(7) 

= 

20.0 

feet 

XCP0(7) 

= 

10.  4 

feet 

forward  of 

transom 

XPO  (8) 

= 

-1.275 

feet 

forward  of 

transom 

XRO  (9) 

= 

-1. 125 

feet 

forward  of 

transom 

XLAERO  (10) 

= 

20.0 

feet 

RSPAN  (9) 

= 

1.  21 

feet 

YSH  (4) 

= 

5.37 

feet 

from  centerline 

AVBMSW  (4) 

= 

0.5 

feet 

WIDTH  (7) 

= 

10.0 

feet 

YPO  (8) 

= 

5.  55 

feet 

from  centerline 

YR(9) 

= 

5.55 

feet 

from  center 

■line 

BEAM  (10) 

= 

10.0 

feet 

ZS(2) 

= 

2.54 

feet 

above  keel 

ZSSI  (5) 

= 

1.875 

feet 

above  keel 

ELMAXS  (5) 

= 

1.875 

feet 

(Program  3 

only) 

ZBSI(6) 

= 

1.  875 

feet 

above  keel 

ELKAXB  (6) 

= 

1.875 

feet 

(Program  3 

only) 

CEHCAB(6) 

= 

1.  875 

feet 

(Program  3 

only) 

XBBW(7) 

= 

10.0 

feet 

BUBHGT  (7) 

= 

1.915 

feet 

ZPO  (8) 

= 

-0.604 

feet 

above  keel 

ZRO  (9) 

= 

-0.208 

feet 

above  keel 

WEIGHT  (2)  =6050. 
THSSI(5)   =   63. 


0    pounds 

4    degrees  (Program  2  only) 
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TKBSI  (6) 
AIXX  (2) 
AIYY  (2) 
AIZZ(2) 
AIXZ(2) 
ALEAK  (5) 
DPSS(5) 
DPBS(6) 
BLEAK  (6) 
RAR2A(9) 

CFSW(4) 

CDSW(U) 

CFSS(5) 

CFBS(6) 

FHCRIT(7)  = 

RASPR(9) 

ETC  (9) 

STHS(16) 

STHP(16) 


=   63.4  degrees  (Program  2  only) 

=2870.0  slug-sguare  feet 

=9320.0  slug-sguare  feet 

=10580.0  slug-sguare  feet 

=-2800.0  slug-sguare  feet 

=   26.0  sguare  feet 

=    1.0  pounds  per  sguare  foot 

=    1.0  pounds  per  sguare  foot 

=    0.1  sguare  feet 


0. 68   sguare  feet 

0.7 

1.  28 
0.  9 
0.9 
0.  556 

2.  15 
0.  167 
0.01 
0.01 


27 


III.  £ITCH_ ANGLE 
(Steady  State  Conditions) 

A.  PROGRAM  2 

1  •  Initial  Runs  and  E val ua t ion 

Using  the  initial  conditions  from  Ref.  1  the  steady 
state  values  of  program  2  were  obtained  as  reported  in  table 
III.  Whereas  the  values  for  draft,  plenum  pressure  and 
thrust  lay  in  the  expected  range,  the  pitch  angle  increased 
with  higher  speed.  This  tendency  was  the  first  point  for  a 
more  thorough  investigation  of  the  changes  made  in  Ref.  2. 


Table_ 

.III 

Steady  State 

Condit. 

Lons 

(Progr 

am  2,  ori 

jinal  version) 

peed 

Pitch 

Draft 

Plenum 

Thrust 

Rudder 

angle 

2: 

cressure 

knots) 

(deg) 

(in) 

(psf) 

(lb) 

(deg) 

10.  0 

1.7^ 

8.  18 

24.4  5 

405.16 

0.0 

12.  5 

1.71 

8.  10 

24.46 

356.42 

0.0 

15.0 

1.74 

7.  97 

24.46 

352.44 

0.0 

17.5 

1  .80 

7.78 

24.46 

373.96 

0.0 

20.0 

1  .90 

7.54 

24.46 

411.18 

0.0 

22.5 

2.01 

7.23 

24.46 

457.84 

0.0 

25.0 

2.15 

6.87 

24.46 

509.57 

0.0 

27.5 

2.29 

6.45 

24.46 

560.98 

0.0 

30.0 

2.47 

6.01 

24.33 

613.63 

0.0 

The  simulation  derived  in  Ref.  1  contained  a  speed 
dependent  correction  term  for  the  x-coordinate  of  the  center 
of  pressure  which  was  assumed  to  move  aft  with  increasing 
speed  and  thus  reducing  the  pitch  angle  [Ref.  6].  This  teem 
had  been  removed  in  Ref.  2  with  the  assumption  that   due   to 
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the  existance  of  a  pressure  wave  [Ref.  7  ]  an  incremental 
wedge  of  volume  has  to  be  added  to  the  plenum  volume  which 
varies  with  speed  and  moves  the  center  of  volume  aft  with 
increasing  speed.  The  according  changes  introduced  in 
subroutines  RHS  and  KAVES  were  expected  to  take  care  of  the 
correct  pitch  angles. 

In  order  to  study  the  effect  of  the  previously 
mentioned  change  -  here  called  water  slope  correction  -  the 
moments  of  the  original  program  2  were  recorded  as  shewn  in 
table  I.  Then,  one  after  the  other,  the  water  slope 
correction  terms  in  the  subroutines  INCON,  BOWSL,  STNSL,  £HS 
and  SIDEWL  were  removed.  Due  to  the  nature  of  these  terms  a 
different  pitch  angle  could  be  expected  only  at  lower  speed. 
In  all  cases,  at  10  knots  for  example,  the  angle  was  lower 
than  the  original  value.  The  biggest  change  in  the 
corresponding  moment  could  be  observed  in  the  bow  seal 
routine  when  the  pitch  angle  dropped  to  1.21  deg. 

The  overall  results  from  these  first  trial  runs 
indicated  that  the  water  slope  correction  provided  a 
neccessary  curvature  in  the  pitch  curve  over  the  speed  range 
of  interest  pointed  out  in  the  following  investigation,  but 
could  not  take  care  of  declining  angles  with  increasing 
speed. 

2-  Selection  of  Pitch  Ancjle  Curve 

As  pointed  out  in  Ref.  5  the  actual  pitch  angles  of 
the  XR-3  are  dependent  on  the  loading  condition  of  the  craft 
and  the  seal  adjustments.  Several  curves  are  presented  from 
which  the  following  two  were  chosen: 

-  The  upper  limit  for   the   pitch   angles   was   set   by   the 
so-called  "ideal"  curve  with  theoretically  derived  values. 

-  The  lower  limit-  was  given  by  the  so-called   "three   turns" 
curve.   These   values   were  actually  measured  on  the  craft 
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after  the  seals  were  adjusted  with  "three  turns"  by   means 

of  a  spindle. 
The  values  for  these  two  curves  and  the  original  pitch 
angles  of  program  2  are  shown  in  table  IV  and  graphed  in 
figure  2.  Additionally  the  values  needed  later  as  initial 
conditions  were  plotted  using  subroutine  COLFIL  after 
modifying  for  thrust  as  a  variable  (plots  1  through  4) . 
From  these  data  the  neccessity  to  correct  the  program  under 
discussion  was  obvious. 

Table  IV 


Soeed 

"Ideal 

(knots) 

(deg) 

10.0 

2. 45 

12.5 

2.  15 

15.0 

1.70 

17.5 

1.25 

20.0 

0.35 

22.5 

0.60 

25.0 

0.30 

27.5 

0.25 

30.0 

0.22 

Steady_S tate_?itch_Anales 
(Original  version) 
"Three  Turns'1 
(deg) 

1.45 
1.  10 
0.80 
0.65 
0.40 
0.30 
0.25 
0.22 
0.20 


Proqram_2 

(deg) 

1. 

74 

1. 

71 

1. 

74 

1. 

80 

1. 

90 

2. 

01 

2, 

15 

2. 

29 

2. 

47 

3.  Changes  and  Justification 

First  Change : 

In  order  to  achieve  the  desired  pitch  angles  one  or 
more  moments  as  listed  in  table  1  had  to  be  changed.  As  it 
was  not  intended  to  modify  the  modeling  of  the  individual 
components  like  seals,  sidewails,  rudder,  propulsion  system 
or  fans,  which  surely  would  endanger  a  qualitative 
comparison   with  program  3,  the  only  moments  left  were  FKAVZ 
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FIGURE  2 


STEADY  STATE  PITCH  ANGLES 


PITCH    ANGLE    (DEGREES) 


2.5    + 


2.0    + 


1.5    + 


1.0    + 


0.5    + 


+----  +  — 

10.0 


15.0 


"IDEAL"    CURVE 
"THREE    TURNS"    CURVE 
SIMULATION    CURVE 
(PROGRAM    2) 

0 


—  +  ,-. 

20.0 


# 

+ 


+  -• 

25.0 


—  +-• 
30.0 


VELOCITY    (KNOTS) 
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and  FMEUE  assuming  the   aerodynamic   forces   found   in   wind 
tunnel  test  being  correct. 

In  figure  3  the  corresponding  forces  and  their  points 
of  attack  are  shown  together  with  the  according  equations  as 
used  in  the  subroutines  INCON  and  RHS.  Following  the 
calculations  it  can  be  seen  that  the  moment  FWAVZ  is 
generated  by  the  bubble  wave  making  drag  FXPWAV,  which 
basicly  is  a  function  of  the  craft's  velocity  and  weight. 
PEBAR  can  be  considered  constant  similar  to  the  pressure 
values  in  table  III.  In  order  to  reach  smaller  pitch  angles 
at  higher  speed  FXPWAV  would  have  to  be  increased 
substantially  with  the  result  that  the  water  slope 
correction  term  wATSLP  and  the  thrust  of  the  craft  attained 
unrealistic  values  at  lower  velocity  where  the  pitch  angle 
tends  to  become  negative  as  well  as  unforseeable  changes  in 
the  dynamic  response  may  be  introduced.  Therefore  this 
possibility  was  rejected. 

As  mentioned  before  the  plenum  pressure  can  be 
considered  constant.  So  the  moment  left  for  a  possible 
pitch  correction,  FMBUB,  could  be  used  in  this  sense  only  by 
changing  its  lever,  XCP,  in  a  speed  dependent  manner.  The 
method  to  determine  the  neccessary  shift  of  the  center  of 
pressure  (C.?.),  in  general,  followed  that  outlined  in  Ref. 
1,  the  resulting  values  are  given  in  table  V  and  graphed  in 
figure  4.  During  this  procedure  it  was  tried  to  match  the 
pitch  angles  of  the  "three  turns"  curve.  It  can  readily  be 
seen  that  there  is  an  almost  linear  relationship  between  the 
location  cf  C.P.  and  the  craft's  velocity.  The  following 
linear  function  used  to  approximate  these  data  points  was 
then  included  in  the  subroutines  RHS  and  WAVES  : 

XCPU  =  XCP  -  0.034*U*0.5921  +  0.3 
This  new  variable  XCPU  was  used   in   the   force   and   moment 
calculations. 
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FIGURE  3 


FORCES  AT  C.G.  AND  C.P 


-  -  XCPO 

XS  -  - ► 
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« 

3 
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I*    f«"---FXPWAv'  '| 


I 


WATERLINE 


KEEL 
(BASEL INE) 


RELATED  CALCULATIONS: 


1.  IN  SUBROUTINE  INCON: 


PWC0N=4.*WEIGHT/(RH0*G*XLBW) 

FNCON=SQRT(G*XLBW) 

XCP=XCPO-XS 
ZCP=ZS-BUBHGT 

U=UO*1.6878 

PBBAR=DELP I 


2.  IN  SUBROUTINE  RHS: 

FN=U/FNCON 
CF=.37/(FN**1. 5655981 ) 

FXPWAV=-PWCDN*PBBAR*CF 
WAT SLP=-FXP WAV/ WEIGHT 

FMBUB=ABPS=MXCP-THETA*Z) 

FWAVZ=-FXPWAV*ZS 
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FIGURE  4 


SHIFT  OF  CENTER  OF  PRESSURE 


(PROGRAM  2) 


XCPO  (FEET) 

I 
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Table  V 


Position,  of  __Cen  ter„pf .  Pressure 
(Program  2,    first  change) 

XCPO 
(feet) 

10.36 

10.25 

10.14 

10.06 

9.95 

9.87 

9.81 

9.75 

9.68 


Speed 

knots) 

10. 

.0 

12 

.5 

15. 

0 

17. 

5 

20. 

0 

22. 

5 

25. 

0 

27. 

5 

30. 

0 

The  updated  steady  state  values  are  contained  in  table 
VI.  Whereas  the  pitch  angles  approached  the  destinated 
values  in  a  satisfactory  way,  the  greater  draft  (and 
conseguently  increased  thrust)  especially  at  higher  speed 
may  be  regarded  a  setback.  But  as  at  this  point  it  could  be 
seen  that  a  similar  correction  would  be  neccessary  for 
program  3  its  impact  on  this  comparison  was  not  considered 
crucial. 

Second  Change^ 

When  calculating  the  moment  due  to  the  bubble  wave 
making  drag,  FHBUB,  program  2  contains  a  correction  term  for 
XCP  in  subroutine  RHS  : 

FHBUB  =  ABPB* (XCP-THfiIA*Z) 
As  Z  is  used  as  a  negative  guantity  (computed  in   subroutine 
IKCON  :   Z=  -ZS+DSO/12.0  ),  FMBUB  would  increase  with  larger 
pitch  angles  provided  A3PB  does  not   vary   much.   Under   the 
small   angle  assumption  for  the  pitch  angle  THETA  this  lever 
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change  generally  is  true,  but  in  the  above  situation  does 
not  represent  the  position  of  C.P.  of  the  XR-3.  Presumably 
the  location  of  C.P.  or.  the  100B  craft  was  almost  underneath 
the  center  of  gravity  (C.G.),  with  ZCP  much  larger  than  XCP 
and  nearly  equal  to  Z.  In  the  case  of  the  XR-3  for  the  given 
data  ZCP  is  only  twice  as  large  as  XCP  and  considerably 
smaller  than  the  average  Z. 

Table  VI 


Speed 

(knots) 

10.  0 
12.5 
15.  0 
17.5 
20.  0 
22.5 
25.  0 
27.5 
30.0 


Steady _ St ate_ Conditions 
(Program  2,  after  first  change) 
Pitch    Draft   Plenum     Thrust 


Rudder 


angle 

pressure 

(deg) 

(in) 

(psf) 

(lb) 

(deg 

1.45 

8.33 

24.34 

405.84 

0.0 

1  .17 

8.  18 

24.46 

358.56 

0.0 

0.94 

8.  10 

24.46 

357.52 

0  .0 

0.76 

7.  99 

24.46 

383.74 

0.0 

0.61 

7.86 

24.46 

428.38 

0.0 

0.48 

7.71 

24.46 

486.60 

0.0 

0.37 

7.54 

24.46 

555.25 

0.0 

0.27 

7.37 

24.46 

632.21 

0.0 

0.19 

7.  18 

24.46 

715.64 

0.0 

Therefore   this   terra  was  removed  and  by  the  following 
function  replaced: 

XCPC  =  SHXYAX(XCPU,  ZCP,  THETA,  PI) 
This   corrected   x-coordinate   was   than  reintroduced  in  the 
subroutines  WAVES  and  RHS,  the  letter  using  the  new  value  to 
calculate  the  moment: 

FMBUB  =  ABPB*XC?C 
The  function  SHXYAX  -  shift  of  x-coordinate  due  to   rotation 
around  the  y-axis  -  is  kept  quite  general  in  order  to  enable 
its  use  for  other  cases.  In  this  subprogram   the   radius   of 
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the   path   that  discribes  the  movement  of  C.P.  is  calculated 
and  its  angle  with  the  X-axis. 

FUNCTION  (X,  Z,  ANGYAX,  PI) 

H  =  SQRT  (X**2+Z**2) 

IF  (X.EQ.0.0)  GO  TO  1 

ARG  =  Z/X 

ANGOLD  =  AT AN  (ARG) 
The   IF-statement   avoids  the  calculation  of  ARG  with  a  zero 
denominator.  Depending  on  the  sign  of  ANGOLD  the  true   angle 
between   H   and   X-axis   is   computed   and  the  (pitch)  angle 
ANGYAX  suttracted: 

IF  (ANGOLD. GE. 0. 0)  GO  TO  2 

ANGNEW  =  ANGOLD  +  PI  -ANGYAX 

GO  TO  3 

1  ANGNEW  =  PI/2.0  -  ANGYAX 
GO  TO  3 

2  ANGNEW  =  ANGOLD  -  ANGYAX 

The    corrected   x-coordinate    is    finally    found    by: 

3  SHXYAX  =  H*COS  (ANGNEW) 

After  introduction  of  this  change  the  simulation 
program  was  re-run  and  the  final  results  recorded  in  table 
VII  and  figure  5.  Though  this  second  change  had  only  little 
influence  on  the  steady  state  conditions  as  measured  it  was 
felt  that  it  should  stay  in  the  simulation  program  because 
of  its  possible  effect  on  the  dynamic  response  where  larger 
pitch  angles  may  occur. 

The  new  moments  are  given  in  table  VIII.  They  were 
attained  under  the  same  run  conditions  near  steady  state  as 
those  of  table  I.  Due  to  the  now  negative  moment  FMBU3  (a 
total  shift  by  about  2000  lb-ft)  and  smaller  pitch  angle  the 
bow  seal  increased  its  moment  considerably-  whereas  the 
sidewalls  reduced  their  effect  as  generally  could  be 
expected . 
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Table  VII 


Steady  .State  Conditions 
(Program  2,  after  second  change) 


Speed 

(knots) 

10.  0 
12.5 
15.  0 
17.5 
20.  0 
22.5 
25.0 
27.5 
30.  0 


Pitch 
ancjie 
(deg) 

1  .50 
1  .21 

0.98 
0.79 
0.63 
0.49 
0.37 
0.27 
0.20 


Draft 

(in) 

8.31 
8.  17 
8.  10 
7.  99 
7.85 
7.71 
7.  54 
7.37 
7.  18 


Plenum 
pressure 
(psf) 

24.36 
24.45 
24.46 
24.  U6 
24.46 
24.46 
24.46 
24.46 
24.46 


Thrust   Rudder 


(lb) 

405.74 
358.42 
357.38 
383.56 
428.20 
486.38 
555.02 
631.98 
715.42 


(deg) 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


Table  VIII 


Summar 

y 

of  Moments 

(Pi 

tch]_ 

(Program 

2, 

.  after  second 

change) 

Component 

Positive 

Negative 

Name 

Moment 
(lb~ in) 

Moment 
(lb-in) 

FMBS 

1227 



FMSS 



81 

FMSW 



2619 

FMRUD 



83 

FMP 

1292 



FMAED 

650 



FMBUB 



109 

FWAVZ 



277 

FMWAV 

0 



Summary 

3169 

3169 
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For  the  further  comparison  with  program  3  the 
variables  pitch  angle,  z-displacement,  plenum  pressure  and 
thrust  (startoard)  versus  time  were  plotted  (plots  5  through 
8)  starting  from  the  initial  conditions  as  given  in  table 
II.  Although  in  this  part  of  the  study  the  steady  state 
conditions  were  of  primary  interest,  some  general 
conclusions  concerning  the  dynamic  behaviour  already  can  be 
drawn.  Plot  5  additionally  indicates  a  more  damped  pitch 
response  due  to  the  first  change. 


Table    IX 

St< 

Bady    State    Condit: 

ions 

(Progi 

:am    3,    original    version) 

S_peed 

Pitch 
ancjle 

Draft 

Plenum 
pressure 

Thrust 

Rudder 

(knots) 

(deg) 

(in) 

(psf) 

(lb) 

(deg) 

10.  0 

0.94 

6.68 

24.88 

521.9 

0.0 

12.  5 

0.68 

6.48 

24.83 

428.0 

0.0 

15.0 

0.47 

6.22 

24.82 

392.8 

0.0 

17.5 

0.28 

5.  94 

24.81 

389.6 

0.0 

20.0 

0.11 

5.65 

24.81 

4  06.3 

0.0 

22.  5 

0.69 

6.23 

24.86 

4  4  9.3 

0.0 

25.0 

0.63 

6.03 

24.  86 

492.7 

0.0 

27.5 

0.56 

5.80 

24.87 

543.6 

0.0 

30.  0 

0.48 

5.46 

24.98 

596.7 

0.0 

Note: 

XS    =     10. 

08    feet 

PROGRAM    3 

1.  Initial  Runs  and  Evaluation 

The  basic  craft  data  used  in  Ref.  3  fortunately  did 
not  differ  from  those  of  Ref.  2  exept  that  the  center  of 
gravity  (C.G)  was   positioned   slightly   more   forward.   The 
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results  of  the  first  simulation  runs  are  summarized  in 
tables  IX  and  X  under  application  of  the  same  initial 
conditions,  the  pitch  angles  are  graphed  in  figure  6.  It  was 
immediately  apparent  that  in  this  program  version,  too,  a 
correction  of  the  pitch  angle  was  neccessary  before  a  first 
comparison  could  be  made. 


Table_X 

Steady  State  Pitch  Angles 

(Program  3,  original  version,  and  Ref.  5) 
Speed    "Ideal"    "Three_Turns^    Program  3 


(knots)   (deg) 


10.0 

2.45 

12.5 

2.  15 

15.0 

1.70 

17.5 

1.25 

20.0 

0.85 

22.5 

0.60 

25.0 

0.30 

27.5 

0.25 

30.0 

0.22 

(deg) 

1.45 

1.10 

0.80 
0.65 
0.40 
0.30 
0.25 
0.22 
0.20 


(deg) 

0.94 
0.68 
0.47 
0.28 
0.  11 
0.69 
0.63 
0.56 
0.48 


2.  Changes  and  Justification 

Zi^st  Change^ 

A  deeper  look  into  program  3  revealed  that  in 
subroutine  RKS  the  speed  dependent  shift  of  C.P.  of  Ref.  1 
was  still  in  use,  whereas  subroutines  INCON,  BOWSL  and  RHS 
already  contained  water  slope  correction  terms  like  Ref.  2. 
As  the  existing  calculation  of  the  shift  of  C.P.  was  not 
satisfactory  (the  cause  for  the  discontinuity  at  22  knots 
remained  unknown),  these  terms  were  removed  from  subroutine 
RHS  and  temporarily  replaced  by  XCP.  Additionally  all   water 
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FIGURE  6 
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slope  correction  terms  were  revised  or  introduced  in 
subroutines  INCCN/  SIDEWL,  RHS,  BOWSL  and  STNSL  as  outlined 
in  Ref.  2  in  all  details.  Besides  these  changes  the  sign 
errors  in  subroutine  PROP  of  Ref.  1  discovered  by  Ref.  2 
were  removed.  A  few  trial  runs,  however,  indicated  only  a 
negligible  effect  on  the  steady  state  values. 

The  results  of  the  following  simulations  are  reflected 
in  table  XI  and  figure  7.  Despite  the  substantial  changes 
made  in  Ref.  3  the  overall  tendency  of  the  craft's  pitch  was 
almost  identical  with  that  of  the  original  version  of  Ref.  2 
and  gave  a  further  proof  that  the  introduction  of  WAISLP 
alone  was  not  sufficient  to  reach  proper  pitch  angles. 

Table_XI 

S  t  e  a d .  y  _  5 1  a  t  e  __  Con  d  i  ti  o  ns 
(Program  3,  after  first  change) 
Speed     Pitch    Draft   Plenum     Thrust   EH^^r 

ancjle         EE^SGure 
(knots)   (deg)    (in)      (psf)     (lb)    (deg) 

10.0  1.97  7.05  24.90  405.67  0.0 

12.5  1.82  7.11  24.92  347.40  0.0 

15.0  1.75  7.09  24.91  334.35  0.0 

17.5  1.75  7.02  24.84  347.26  0.0 

20.0  1.78  6.93  24.71  377.27  0.0 

22.5  1.85  6.76  24.61  415.74  0.0 

25.0  1.97  6.49  24.51  459.64  0.0 

27.5  2.12  6.13  24.45  504.41  0.0 

30.0  2.15  5.93  24.35  555.91  0.0 

Note:  XS  =  10.08  feet 

Furthermore  the  center  of  gravity  was  re-positioned  in 
order  to  match  the  conditions  used  in  Ref.  2  (see  II. C)  . 
This   seemed  permissable  as  both  Refs.  2  and  3  used  the  same 
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FIGURE  7 
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craft  configuration  and  attained  a  shift  of  C.G.  by  adding 
some  load  to  the  craft.  The  new  location  mainly  effected  the 
pitch  angle  as  indicated  in  table  XII.  The  simulation  curve 
of  figure  7  was  lifted  up  by  an  almost  constant  amount,  but 
the  general  shape  was  not  influenced.  Therefore,  the  task 
remained  to  introduce  an  additonal  correction  term. 

Table  XII 


iii£MY_State_Conditions 
(Program  3,  after  first  change  and  shift  of 


G.) 


Speed 

(knots) 

10.  0 
15.  0 
20.0 
25.  0 
30.0 


Pitch 
ancjle 
(deg) 

2.  13 

1  .90 
1  .92 
2.  11 
2.22 


Draft 

(in) 

7.  03 
7.07 
6.  97 
6.46 
6.  10 


Plenum 
pressure 
(psf) 

24.87 
24.87 
24.64 
24.43 
24.31 


Thrust   Rudder 


(lb) 

404.73 
332.99 
375.79 
455.68 
558.27 


(deg) 

0.0 
0.0 
0.0 
0.0 
0.0 


At  this  point  the  corresponding  moments  near  steady 
state  were  recorded  in  table  XIII  which  enables  a  direct 
comparison  with  program  2  under  similar  conditions  (see 
table  I)  .  The  reduced  draft  and  plenum  pressure  as  well  as 
the  changes  made  in  Ref.  3  can  be  considered  to  account  for 
the  different  moments  FMBS,  FMSS,  FMSW  and  FMBOB.  Again  the 
same  variables  as  in  program  2  were  output  (plots  9  through 
12)  . 
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Table_XIII 

Su  m  mar  x._of  foments (pitch) 

(Program  3,  original  version) 
Component     Positive       Negative 
Name  Mo^HSBi  Moment 

(lb-in)         (lb-in) 

FMBS  850             

FMSS            324 

FMSW            3694 

FMRUD           83 

FMP  1049             

FMAED  651              

FMBUB  1831             

FWAVZ           280 

FMWAV  0             

Summary        4381  4381 

Second  Change^. 

The  following  modification  was  carried  out  the  same  way 
as  the  first  one  of  program  2  (see  III. A. 3)  .  The  resulting 
values  are  represented  in  table  XIV  and  graphed  in  figure  8. 
As  a  linear  curve  fit  did  not  suffice  a  satisfactory 
approximation  any  longer  a  guadratic  function  was  developed 
with  final  form: 

XCPU=XCP+0.001975* (U* 0.5921-30.0)  ** 2- 0.9  74 
This   new   variable  was  used  the  same  way  to  calculate  force 
and  moment.  The  resulting  pitch   angles   besides   the   other 
steady  state  values  are  shown  in  table  XV. 
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FIGURE    3 
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Table, 

_XIV 

Position   of 

Cen  t 

?r 

of    Pressure 

(Program 

3,    secor 

id    change) 

S£>eed 

XCPO 

(knots) 

(feet) 

10.0 

10.22 

12.5 

10.05 

15.0 

9.85 

17.5 

9.71 

20.0 

9.56 

22.5 

9.50 

25.0 

9.46 

27.5 

9.42 

30.0 

9.42 

Table  XV 

Steady  S t a t e _ Con d i t i o ns 
(Program  3,  after  second  change) 

S£eed     Pitch    Draft   Plenum     Thrust  Rudder 

(knots)   (deg)    (in)      (psf)     (lb)  (deg) 

10.0      1.36     7.11     24.89     407.85  0.0 

12.5      1.03     6.89     2a. 86     348.84  0.0 

15.0      0.78     6.61     24.86     334.54  0.0 

17.5      0.59     6.31     24.85     345.26  0.0 

20.0      0.44     6.05     24.84     372.28  0.0 

22.5      0.33     5.81     24.84     410.41  0.0 

25.0      0.26     5.61     24.84     457.48  0.0 

27.5      0.23     5.44     24.85     511.94  0.0 

30.0      0.20     5.17     24.97     568.62  0.0 
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liii£d  Chanc]iej_ 

The  last  alteration  reflected  the  shift  of  C.P.  due  to 
the  craft's  motion  and  was  done  by  introducing  -che 
corresponding  function  SHXYAX  which  is  described  in  detail 
in  III. A. 3.  table  XVI  and  figure  9  show  the  final  results 
with  pitch  angles  slightly  off  the  anticipated  values.  As 
this  was  only  marginal  a  new  curve  fit  was  not  done,  but 
once  knowing  the  quadratic  relationship  it  would  De  a 
relatively  short  and  straight  forward  procedure  (see  V). 


Table_ 

_XVI 

Steady  Sta. 

te 

Conditd 

.ons 

(Progra 

m  3,  af 

ter  third 

change) 

peed 

Pitch 

Draft 

Plenum 

Thrust 

Rudder 

angle 

21 

cessure 

knots) 

(deg) 

(in) 

(psf) 

(lb) 

(deg) 

10.  0 

1.38 

7.  16 

2  4.84 

407.95 

0.0 

12.5 

1  .05 

6.  93 

24.84 

384.90 

0.0 

15.0 

0.80 

6.64 

24.84 

334.64 

0.0 

17.5 

0.60 

6.  33 

24.84 

345.32 

0.0 

20.  0 

0.45 

6.  06 

24.3  4 

372.36 

0.0 

22.5 

0.34 

5.82 

24.84 

410.52 

0.0 

25.0 

0.28 

5.  62 

24.84 

457.54 

0.0 

27.5 

0.24 

5.45 

24.84 

5  12.04 

0.0 

30.  0 

0.2  4 

5.30 

24.84 

568.60 

0.0 

Some  of  the  new  steady  state  moments  have  changed 
considerably  as  can  be  seen  from  table  XVII.  Due  to  the 
further  aft  moving  C.P.  most  of  the  compensating  moment  is 
picked  up  by  the  bow  seal.  It  must  be  notified  that,  in 
general,  the  total  moments  of  the  corrected  program  3  were 
significantly   larger   than   those   of    program    2    after 
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FIGURE  9 
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modification   and  would  have  had  greater  o   ect  on  the  shear 

forces  and  bending  moments  v?hich  are  not  considered  in   this 

simulation (see   V).   The   listed  variables  were  also  plotted 

for  comparison  studies  (plots  13  through  16). 

Table_XVII 

Su  m  ma  r_Y_of_Mo  merits (£itch]_ 

(Program  3,    after  third  change) 
Comp_onent     Positive       M^£aiiv§ 
Name         H2!ii:§nt         Moment 
(lb-in)         (lb-in) 

FMBS  2985             

FMSS  ---  257 

FMSS            2060 

FMRUD           83 

FMP  1138             

FMAED  651             

FMBUB           2094 

FWAVZ           280 

FMWAV  0              

Summary        4774  4774 

C.  COMPARISIOK  AND  RESULTS 

Although  the  data  so  far  obtained  allowed  some  comparison 
and  conclusions  additional  information  was  neccessary  to 
demonstrate  the  basic  differences  between  the  two  program 
versions.  In  table  XVIII  the  interesting  variables  are 
listed  which  were  obtained  in  typical  runs  at  20  knots  after 
all  changes  mentioned  had  been  introduced.  Also  the  time 
histories  for  pitch  angle,  draft,  plenum  pressure  and  thrust 
(starboard)  give  a  good  basis  for  a  qualitative   evaluation. 

Starting  from  these  trajectories  a  significant  difference 
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can  be  seen  in  the  pitch  angle  oscillations  and  the  damping 
factor  which  will  be  investigated  more  thoroughly  in  the 
following  chapters.  The  higher  thrust  in  program  2  is 
certainly  due  to  the  greater  draft  so  that  under  these 
consideration  both  programs  reached  reasonable  steady  state 
conditions. 

Table_XVIII 

§teady_Sta te_?er forma nce_ Da ta 
(Modified  programs  2  and  3) 

Variables  Mni^§  P.L21'£§.Q-Z  HLQ.9lIL$.E-1 

Speed  kt  20.0  20.0 

Pitch  angle  deg  0.63  0.45 

Time  constant  sec  3.36  3.06 
Freguency  of 

oscillation  Hz  0.48  0.89 

Draft  in  7.85  6.06 

Thrust  lb  428.20  372.36 
Air  flow  rates 

Bow  fans  cu  ft/sec  125.3  20.3 

Main  fans  cu  ft/sec  135.2  22.0 

Stern  fans  cu  ft/sec  125.3  12.9 

total  cu  ft/sec  385.8  55.2 

Fan  power  required     hp  17.61  2.55 

Actual  fan  power  hp  ■  8.93 

Plenum  pressure  psf  24.46  24.84 

Leakage  flow  rates 

Bow  seal  cu  ft/sec  12.9  10.4 

Sidewalls  cu  ft/sec  0.0  0.0 

Stern  seal  cu  ft/sec  372.9  44.8 

total  cu  ft/sec  385.8  55.2 

Plenum  volume  cu  ft  288.17  318.75 

Nevertheless,  program  3   represents   more   realistic   air 
flow   and   leakage   rates   in   connection  with  the  available 
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plenum  volume  as  indicated  in  table  XVIII  (The  difference  in 
the  volujiina  basicly  a  consequence  of  different  draft  and 
pitch  angles,  too)  .  \-'::e  changes  introduced  had  negligible 
effect  on  the  computation  time  in  both  cases,  but  program  3 
needed  slightly  less  CPU  time  due  to  its  better  damping 
factor. 
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IV.  ROLL_ANGLS 
(Dynamic  Response) 

A.  INITIAL  RUNS  AND  EVALUATION 

In  order  to  study  the  dynamic  behaviour  of  the  simulation 
models  several  possibilities  were  given  to  introduce  the 
neccessary  disturbance.  Regular  or  irregular  sea  states  in 
tabular  form  or  by  superimposed  sine  waves  could  have  been 
entered  using  subroutine  WAVES.  Another  way  was  to  change 
speed  by  means  of  a  thrust  map  of  data  points  converted  into 
a  time  history  in  subroutine  PROP.  Finally  rudder  angles 
could  have  been  applied  following  a  straight  run  or  by 
varying  angles  given  as  data  points,  too.  The  choice  was 
determined  primarily  by  the  idea  that  besides  the 
qualitative  comparison  of  mathematical  models  a  future 
validation  of  the  results  obtained  should  be  made  possible. 
Therefore,  an  easy-to-measure  calm  water  situation  was 
prefered.  As  additionally  new  changes  effecting  the  roll 
damping  were  to  be  incorporated  the  excitation  through  the 
speed  did  not  seem  very  promising.  Once  the  steady  state 
conditions  being  found  and  used  for  initialisation  a  step 
rudder  input  (or  a  steep  ramp  input  to  come  closer  to  the 
real  case)  seemed  the  best  solution  offering  a  fine  control 
of  the  disturbance  magnitude. 

At  a  speed  of  20  knots  -actually  the  corresponding 
initial  thrust  was  kept  constant  -  a  rudder  angle  of  35 
degrees  was  found  to  give  easily  observable  pitch  and  roil 
angle  excursions  for  both  models  whereas  larger  values 
eventually  caused  water  contact  with  the  top  of  the  plenum 
chamber  and  stopped  the  program  execution  (controlled  by 
subroutine  SIDEWL) . 

The  first  simulated  turns  under  these  conditions  were 
recorded   in   plots   17   through   24    for    both    programs 


54 


respectively  (pitch  and  roll  rates  were  changed  into  output 
variables  in  subroutine  COLFIL  and  graphed  in  order  to  ease 
the  estimation  of  damping  factor  ratios) .  As  previously 
learned  program  2  showed  a  less  damped  pitch  angle,  but  its 
roll  angle  response  -  after  removal  of  the  superimposed  low 
frequency  component  of  undetermined  origin  -  is  slightly 
more  damped  than  that  of  program  3.  The  speed  dependent 
shift  of  XCPO  therein  may  also  be  the  cause  for  an  increased 
tendency  towards  self  contained  pitch  angle  oscillation. 

B.  ORGANIZATION  0?  SUBROUTINE  SIDEWL 

Before  the  next  changes  are  explained  a  short 
introduction  into  the  subroutine  program  SIDEWL  is  presented 
in  order  to  give  a  better  understanding  of  the  possibilities 
and  limitations  of  these  alterations. 

The   hydrodynamic   and   hydrostatic   forces   and   moments 

acting   on  the  sidewalls  are  found  by  application  of  slender 

body  theory.  This  theory  is  based  on  motion  of  a  foil   in   a 

viscous   medium.   The   methods   by   which   the  lift  and  drag 

forces  associated  with  moving  underwater   bodies   are   found 

are   primarily   linear   in  nature  (Ref.  *J  and  8),  but  heave, 

pitch,   and   roll   motions   cause   the    most    significiant 

nonlinear    variations.    Additionally    there   are   certain 

nonlinear  terms  due  to  cross-flow  drag,  which  are   important 

for  the  case  of  very  low  aspect  ratio  lifting  surfaces  as  in 

the  case  of  the  XR-3  sidewalls.   The   sidewall   forces   will 

then   be   composed   of   terms   due   to  sidewall  buoyancy  and 

slender   body   hydrodynamic   reactions   and   the   effect   of 

cross-flow  drag: 

F  =  F  +  F 

sidewalls     slender  body     cross-flow  drag 

In  subroutine  SIDEWL  first  the  draft  of  or  the  gap 
beneath  the  sidewalls  is  calculated.  The  draft/gap  is   found 
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for  each  X-direction  section  after  it  has  been  corrected  for 
craft  roll  and  pitch  angles  and  waves  including  water  slope 
effect.  In  the  presence  of  gaps  the  total  sidewall  leakage 
area  ALSK  is  computed  by  summing  the  products  of  gap  heights 
and  sidewall  section  lengths  for  both  sidewalls. 

In  the  following  force  and  moment  calculation  the 
cross-flow  drag  terms  and  slender  body  theory  terms  are 
computed  in  separate  sections  and  finally  summed  to  give  the 
total  values.  In  the  case  of  the  roll  moment  FK  this 
calculation  is  of  the  form: 

FK  =  (FZH(2)-FZH (1) ) *YS W-ZS*FY+FKD 
The   slender  body  theory  terms  FZH  for  each  sidewall  (1-port 
side,  2=starboard)  are  given  by: 

FZH  =  -G*BC0-U*U*A33S*THETA-U*A33S*W+Q*U 
*  (A3  3S*XSS-3C2)-U*A3  3S*P*YLSW 
where  the  expressions  used  are  defined  as  follows: 
BCO   =  mass  of  displaced  water  by  the  jth  sidewall  element 
A33S  =  vertical  added  mass  of  jth  sidewall  element 
XSS   =  -XS 

BC2   =  total  sidewall  vertical  added  mass 
YLSW  =  lateral  distance  of  jth  sidewall  element  from  the 

centerline 
These  components  are  either  calculated  in  the  subroutine 
itself  or  looked  up  in  the  sidewall  tables  computed  in 
subroutine  SIDETAB.  Further  details  concerning  the  set-up 
of  these  tables  are  described  in  Ref.  4,  the  resulting 
values  are  stored  in  COMMON-block  KAVTAB.  The  total  lateral 
force  FY  is  defined  by: 

FY  =  FYH  (1)  +FYK (2)  -FYD 
It  represents   the   hydrody naniic   lateral   forces   plus   the 
corresponding  cross-flow  drag  component  derived  as  follows: 

FYH  =  -A22S*U*(V+XSS*R-ZS*P) 
where  A22S  is  the  lateral  added  mass  at  the  stern  given  by: 

A22S  =  (RHO*.4*PI*DSS**2)/2 
with   DSS  beeing  the  computed  sidewall  draft.  The  final  term 
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FKD  is  entirely  due  to  cross-flow   drag   on   both   sidewalls 
each  of  them  defined  by: 

FKD  =  FYD*  (ZS-DSWAV/2) 
where  FYD  is  the  lateral  force  due  to  cross-flow  drag  and 
DSWAV  is  the  draft  of  the  jth  sidewall  section  corrected  for 
waves.  In  the  computation  of  FK  the  lateral  force  FY  is 
used  without  the  drag  component  which  is  added  later.  In 
general  it  must  be  said  that  all  forces  more  or  less 
influence  each  other  so  that  a  change  in  the  roll 
characteristics  neccessarily  must  alter  the  pitch  and  yaw 
behaviour  or  vice  versa  as  could  be  seen  when  intrducing  the 
pitch  angle  corrections  or  in  the  first  turn  runs  where  the 
pitch  angles  increased  significantly. 

C.  CHANGES  AND  JUSTIFICATION 

The  following  changes  were  recommended  by  Bentson  [Ref. 
9].  They  consisted  of  two  major  blocks  and  a  minor  addition 
carried  out  in  subroutine  SIDEWL.  The  letter  dealt  with  the 
vertical  added  mass  at  the  stern,  A22S,  which  is  recomputed 
for  the  case  of  negative  pitch  angles  using  the  corrected 
draft  at  the  bow,  DRBOW: 

DRBOW=DSS~ (XX (J,N  +  1) -XSS) *THETA 

IF  (DRBOW.LT.0.0)  DRBOH=0.0 

A22S=  (RKO*.8*PI*DSS*#2)  /2 . 

IF  (THETA.LT.O. 0)  A22S= . 8*RHO*PI*DRBOW*DRBOW/2. 
XX  is  the  longitudinal  distance  of  the  nth  sidewall  station 
from  the  bow.  The  factor  .8  in  the  calculation  of  A22S  was 
changed  to  .4  in  order  to  remain  consistent  with  the  model 
given  based  on  the  definitions  derived  in  Ref.  4.  The 
effect  of  this  change,  however,  was  not  observable  because 
the  simulated  pitch  angles  were  most  of  the  time  positive 
except  a  few  small  negative  transient  peaks. 

The  next  major  change  was  a  correction  of  the  vertical 
force  FZH.   As  the  outside  edge  of  the  submerged  portion   of 
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the  sidewall  is  not  vertical  near  the  bow  an  additional 
upright  force  component  is  generated  due  to  the  deadrise 
projection  of  the  lateral  force  FYH,  in  algebraic  terms 
defined  by  the  product  of  the  lateral  force  and  the 
cotangent  of  the  deadrise  angle  (average  angle  between  the 
horizontal  and  the  inclined  outside  edge  of  the  mentioned 
sidewall  portion) .  In  FORTRAN  IV  language  the  computation 
proceeds  as  follows: 

CTNDR=0. 0 

IF  (DSS.LE.0.0)  GO  TO  22 

CTNDR= (BS-BB(1) ) /DSS 

IF  (THETA.LT.0.0)  CTNDR=0. 39391 
22  CONTINUE 

FZK  (J)  =FZH(J)  +Pn1*FYH  (J)  *CTNDR 
BS   and   3B(1)   are   the   lower  and  upper  beam  widtas  of  the 
submerged  sidewall  portion,  CTNDR  stands  for  cotangent,   and 
PW1   is   a  unity  sign  factor  depending  on  the  sidewall  under 
discussion. 

The  last  addition  had  the  purpose  to  improve  the  roil 
damping  of  the  model.  It  was  developed  from  the  idea  that 
when  the  craft  is  rolling  the  sidewails  generate  asymmetric 
waves  depending  also  on  speed  which  take  energy  out  of  the 
ship  and,  therefore,  have  a  damping  effect  on  the  roll 
motion.  Without  discussing  each  single  step  of  this  rather 
involved  change  the  general  procedure  was  to  evaluate  new 
draft  numbers  corrected  for  these  waves  and  from  there  to 
find  correction  terms  for  the  vertical  added  sidewall  mass 
BC2.  The  previously  determined  roll  moment  FK  was  finally 
adjusted  by  subtracting  the  resulting  counter  moments 
including  an  experimental  factor.  The  complete  change  was 
done  as  follows: 

DSS=Z+ZS-XSS*THETA 

Z0R1= (SIGN (1. , DSS) +1.) /2. 

DS£=DSS*Z0R1 

DS=Z+ZS 
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DSE=DS- (XREF-XS)  *THETA 

ID=1.+  (DSR+12.-SDS)  /DDS 

ID=MAXO  (MINO(ID,NDS) ,1) 

DDSR=  (ID-1)  *DDS+SDS 

ID1  =  MIN0  (ID+1 rNDS) 

DID= (DSH*12.-DDSR)  /DDS 

BC2=AC2(1,ID,IP) 

BC2=BC2+DID* (AC2 ( 1 , ID1 , IP) -BC2)  +  DIP* (AC2 (1,ID,IP1) 
-BC2  +  DID*  (AC 2  (1,10  1,1? 1) -AC 2  (1,ID,IP1) 
-AC2  (1,ID1,IP)  +  BC2) ) 

FK=FK-16.*YSW*YS*J*BC2*P/PI 
For  the  shear  and  bending  moment  calculations  not   contained 
in  this  program  version  the  vertical  force  components  of  the 
individual  sidewalls  were  also  corrected: 

FZI]  (1)  =  FZH  (1)  +8.*YSW*BC2*P/PI 

F2H  (2)  =FZH  (2)  -8.  *YSW*BC2*P/PI 

Reference  9  finally  indicated  the  option  to  include  the 
calculation  of  the  spray  drag  of  the  sidewalls  when 
operating  at  high  speed  with  minimum  draft,  whereas  this 
function  originally  was  developed  for  the  100B  craft  with 
considerably  smaller  draft  values  down  to  one  inch  at  about 
80  knots  the  XR-3  operates  with  five  to  six  inches  or  more 
and,  therefore,  the  effect  of  the  spray  drag  was  considered 
negligible. 

D.  COMPARISON  AND  RESULTS 

Starting  with  program  2  the  new  roll  and  pitch 
characteristica  were  recorded  in  plots  25  through  28. 
Besides  a  reduced  overshoot  of  the  pitch  angle  by  about  20 
per  cent  the  roll  angles  reached  only  about  half  of  their 
original  peak  values  while  the  non-linearities  of  tne 
sidewall  changes  b.ecame  quite  apparent.  Therefore,  the  time 
constants  could  not  but  be  estimated  and  seemed  to  have 
increased  slightly  for  both  motions.   Although   the   overall 
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response  could  be  considered  an  improvement  "an  order  of 
magnitude  difference  in  the  roll  damping  value"  as  pointed 
out  in  Pef.  9  and  meant  for  the  Bell  100-ton  test  craft 
could  not  be  experienced  for  this  model. 

For  program  3,  whose  time  histories  were  graphed  in  plots 
29  through  32  the  changes  introduced  were  not  as  well 
favourable.  The  roll  angle  overshoot  reduced  only  by  20  per 
cent  whereas  an  equally  diminshed  pitch  angle  seemed 
considerably  less  damped  and  showed  the  increasing  trend  for 
self-contained  oscillation  as  indicated  earlier.  However,  it 
must  be  noticed  here  that  program  3  originally  was  not  quite 
as  underdamped  as  program  2  and,  therefore,  could  not  be 
effected  by  these  changes  by  the  same  amount. 

In  order  to  attain  an  optimal  model  based  on  the  previous 
results  and  considerations  the  two  major  changes  made  in 
subroutine  SIDEWL  and  their  individual  impact  on  the 
simulated  performance  were  studied  more  thoroughly  after 
introducing  a  few  modifications.  In  the  calculation  of  the 
vertical  force  due  to  deadrise  projection  of  FYH  the 
statement 

FZH  (J)=FZH  (J)  +PM1*FYH  (J)  *CTNDR 
was  replaced  by  the  following  expressions: 

FZHOLD{J)=FZH  (J) 

FZHDRP(J) =PM1*FYH (J) *CTNDR*PR0M01 

FZH  (J)  =FZH  (J)  +FZHDRP  (J) 
These  statements  allowed  to  print  the  individual  force 
components  if  desired  and  to  vary  the  component  value  due  to 
deadrise  projection  by  means  of  the  program  card  deck  using 
one  of  the  modification  settings  (see  Appendix  A)  .  Similary 
the  roll  damping  computation  was  re-organized: 

FK0LD=FK 

FK=FK-PROM2*YSW*YSW*BC2*P/PI 

FZH(1)  =FZH(1)  +  PRO M 02/2.  *YSW*BC2*P/PI 

FZH (2) =FZH (2)  -PR0M02/2 . *YSW*BC2*P/PI 
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Finally  a  WRITE-statement  was  added  using  PR0M03  as  print 
switch  for  VAL(1j,  FZH0LD(1),  FZHOLD  (2)  ,  FZHDRP(1), 
FZHDRP(2),  FKOLDr  and  FK.  The  influence  of  various  settings 
of  PR0M01  and  PR0M02  -  here  simply  called  deadrise  factor 
and  roll  damping  factor  -  could  be  studied  on  the  following 
plots  33  through  60,  the  meaning  of  the  R.D.  number  is 
summarized  in  table  XIX. 
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Using  program  3  in  these  simulation  runs  it  could  be 
shown  that  the  response  was  improved  by  lowering  the 
deadrise  factor  and  increasing  the  roll  damping  factor  as 
graphed  in  plots  57  through  60  under  R.D.  number  8.  These 
values  were  chosen  to  indicate  an  approach  method  only,  more 
extended  simulations  with  a  finer  factor  spacing  could 
possibly  result  in  an  even  better  damped  behaviour.  This 
step,  however,  should  be  done  in  conjunction  with  a 
validation  of  this  model  in  comparison  with  measured  test 
data. 
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v  •  CONC  LUS  IONS_ ft  ND_RECOil  H  EN  D  ATION S 

A  great  deal  of  time  -  more  than  anticipated  at  the 
onset  -  was  required  for  familiarization  with  the  computer 
models  and  for  the  implementation  of  comprehensive  and 
reasonable  initial  trim  conditions.  For  future  studies 
Appendix  A  and  the  example  program  printout  as  well  as  the 
presentation  of  many  tabulated  data  may  help  to  reduce  this 
introductory  phase  effectively. 

The  simulated  roll  and  pitch  motions,  in  general,  were 
improved  although  some  of  the  changes  made  need  further 
investigation  in  connection  with  validation  tests  which 
should  be  performed  with  program  3  as  the  presently  more 
advanced  version.  In  case  of  future  changes  effecting 
longitudinal  and  athwartships  stability  it  is  recommended 
that  the  shift  of  center  of  pressure  is  verified  as  outlined 
in  Appendix  A  before  concluding  tests  are  undertaken. 

These  further  studies  may  concentrate  on  the 
investigation  of  the  moment  due  to  the  product  inertia  term 
AIXZ  which  has  some  importance  in  maintaining  stability  in 
roll  and  pitch  for  larger  rudder  angles.  The  use  of  the 
discrete  mass  distribution  option  (Block  2)  can  be  helpful 
for  this  case. 

Anotner  field  of  interest  can  be  the  air  pressure 
distribution  in  plenum  chamber  and  both  seals  and  its 
influence  on  transient  motions. 

Perhaps  the  most  significant  improvements  may  arise  from 
a  comprehensive  study  of  the  modelling  of  the  sidewalls. 
Besides  the  calculation  of  the  vertical  force  due  to 
deadrise  prcjectioh  of  the  lateral  force  -  the  sidewall 
wetted   surface   shape   of   the  Bell  100-ton  craft  is  almost 
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uniform  from  tow  to  stern,  whereas  the  XR-3  is  more  formed 
like  a  wedge  widening  at  the  stern  -  the  impact  of  fiat 
portions  on  the  bottom  of  the  sidewalls  as  indicated  in  Ref. 
3  should  be  investigated  as  they  become  effective  planing 
surfaces  at  higher  speed  (20  -  25  knots)  but  are  not 
considered  in  the  simulation  program  so  far.  The  need  to 
lower  the  computed  cotangent  value  indicates  that  the 
present  calculation  of  the  deadrise  projection  does  not 
enough  consider  these  structural  differences.  A  pitch  down 
due  to  a  more  negative  sidewall  moment  would  also  reduce  tne 
neccessary  shift  of  the  center  of  pressure  and  thus  the 
destabilizing  effect  on  the  transient  response. 
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APPENDIX       A 


SES    MOTIONS    AND    LOADS    PROGRAM 


USERS    MANUAL 


(short    form) 
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INTRODUCTION 


This  users  manual  is  a  short  form  of  that  given  in 
reference  (1)  which  should  be  referred  to  in  case  of  severe 
programming  difficulty  and  for  a  deeper  study  of  tne 
subroutine  programs  involved.  The  version  given  here 
reflects  the  changes  made  in  the  meantime  by  Refs.  2  and  3 
and  by  the  author  of  the  thesis  study  the  manual  is  attached 
to.  Its  objective  is  to  introduce  the  reader  to  the  general 
program  flow,  to  the  neccessary  requirements  to  establish  an 
own  input  deck  for  various  types  of  craft  runs  and  to  the 
output  options. 
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PROGRAM  ORGANISATION 

Within  a  given  run  the  program  proceeds  as  follows   (see 
Figure  10) : 

-  First  the  main  program  calls  the  input  and  initialization 
subroutine  INCON.  The  program  then  proceeds  to  calculate 
the  motion  time  histories  by  calling  the  various  modular 
subroutines  and  the  integration  routine  for  as  many  time 
steps  as  is  neccessary  to  complete  the  run.  During  the 
course  of  the  calculation  the  various  subroutines  will 
output  those  variables  selected  by  the  print  option 
switches  at  the  appropriate  input  print  interval.  In 
addition,  data  Jiay  be  written  on  scratch  files  to  be  used 
for  summarizing  the  output  and  doing  bending  moments  and 
shear  calculations  (at  the  present  time  the  shears  and 
moments  subroutine  has  been  removed) . 

-  When  the  run  is  completed,  the  program,  depending  on  the 
options  chosen,  may  then  print  and/or  plot  the  output 
summary.  After  these  steps  are  completed  the  main  program 
returns  to  INCON  to  read  data  for  the  next  case. 

-  The  above  procedure  is  repeated  until  all  cases  are 
finished. 
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FIGURE  10 


PROGRAM  FLOW  DIAGRAM 
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PROGRAM  REQUIREMENTS 


As  run  on  an  IBM  360  computer  under  OS/360MVT  version 
(Release  20.6) ,  on  W.  R.  Church  Facility,  NPS,  the  program 
requirements  are  as  follows: 

TOTAL  CORE 
Program  length  201,456 


10 


Core  size  reguireed  for  execution   250,000 


10 


The  total  core  required  for  a  given  step  includes  the 
program,  system  routines,  data  storage,  etc.  The  time 
required  to  run  a  single  case  depends  on  the  nature  of  the 
case  and  the  length  of  time  to  be  simulated.  The  ratio  of 
required  execution  time  to  real  time  has  been  found  to  be  as 
low  as  1/3  for  some  calm  water  cases.  Higher  ratios  occur 
where  the  transient  curves  are  rather  sharp  as  during  the 
turn  maneuvers  in  this  study  where  20  seconds  run  time 
needed  about  six  minutes  using  FORTRAN  G  compiler.  However, 
Mitchell  [  Ref.  10]  found  an  average  improvement  of 
thirty-seven  per  cent  CPU  time  by  compilation  in  FORTRAN   H. 

The   program   currently   uses  the  following  FORTRAN  Data 
Set  Reference  Numbers  (Unit  Numbers) : 
Unii  Number     Usacje 

1  Scratch  file  for  COLFIL  plotting  package. 

2  Scratch  file  for  lateral  plane  output 

3  Summary 

5  Input 

6  Output 

10-15   Sidewall  Integral  Taole  data 
IBMF   Intermediate  file  containing  time  histories 
required  by  bending  moment  calculations.  The 
unit  number  IBMF  is  read  in  using  the  01501 
input  card  not  contained  in  the  present 
program  version. 


69 


INPUT 


The  input  logic  is  subdivided  into  twenty-two  "blocks". 
Each  block  being  for  a  particular  subroutine  or  function 
e.g.  stern  seal,  sidewall,  etc.  Switching  to  the  proper 
input  block  is  done  by  means  of  a  control  tag  on  the  input 
cards.  Once  control  has  been  transferred  to  a  particular 
block  additional  input  may  be  read  by  specific  format 
statements  if  necessary,  however,  the  card  which  contains 
the  control  tag  may  also  be  used  to  input  up  to  seven 
variables  for  the  selected  block.  The  format  of  the  input 
card  used  for  switching  is  described  later,  a  sample  data 
deck  for  the  last  simulation  run  (R.D.  number  8)  is  copied 
at  the  end  of  this  appendix. 

The  control  tag  used  to  transfer  control  tc  the 
appropriate  block  is  of  the  form  OMN  where  M  and  N  are 
integer  constants,  i.e.,  if  the  control  tag  is  0  15  the 
program  will  go  to  statement  1500. 

The  option  tag  is  loaded  into  a  variable  called  I0PT  and 

allows  different  branches  of  logic   within   a   given   input 

block.   If   no   branches  are  used  within  the  input  block  the 
option  tag  may  be  zero. 

The  input  variables  are  stored  in  an  array  called  TEMP 
and  may  be  set  to  a  given  parameter  by  using  an  arithmetic 
replacement  statement  in  the  appropriate  block  (see  various 
blocks  in  attached  listing  for  examples) .  Print  switches 
control  the  printing  by  each  subroutine  of  its  forces  and 
moments.  This  option  is  very  useful  in  debugging  the 
program.  Normally  it  provides  too  much  output  and  is  too 
cumbersome. 
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OUTPUT 

The  output  generally  consists  of  the  following  parts: 

-  The  data  input  deck  is  first  copied  headed  by  the  title 
"LISTING  OF  INPUT  DECK".  Each  block  and  option  is  printed 
as  actually  punched  onto  the  cards.  These  data  are  written 
in  the  main  program  prior  to  any  test  and,  therefore, 
appear  in  all  runs. 

-  After  the  data  are  read  by  subroutine  INCON  a  more 
explanatory  summary  is  output  entitled  "SES  MOTIONS  AND 
LOADS  PROGRAM"  followed  by  the  description  written  in 
block  18.  This  table  contains  most  of  the  input  data  and 
first  computation  results  like  inertia  matrix  and  xx-  and 
yy-arrays  of  sidewalls  and  seals.  In  cases  of  input  errors 
a  corresponding  message  follows  and  execution  stops. 
Likewise  all  other  problems  preventing  a  proper  completion 
of  the  current  run  are  indicated  at  this  place. 

-  If  plots  are  requested  they  are  either  printed  in  the 
following  part  or,  for  CALCOMP  plots,  their  titles  and 
completion  are  summarized  using  subroutine  COLFIL. 

-  In  the  final  portion  the  tabulated  variables  as  desired 
are  summarized  in  one  or  two  divisons  followed  by  the 
message: 

COMPLETED  ALL  RUNS 

On  the  following  page  the  possible  variables  for  these 
plots  or  summaries  are  listed.  A  change  to  other  variables 
needs  only  moderate  program  modifications  in  subroutines  RHS 
and  COLFIL,  eventually  in  some  others,  too. 
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Listing  of.  Output  Variables 


The  names  mentioned  below  are  used  in  subroutine  COLFIL, 
but  may  have  different  names  in  other  subroutine  programs. 

Description 

Independent  variable  time,  seconds. 

Wave  height,  feet. 

Z  displacement,  inches. 

Pitch  angle,  degrees. 

Plenum  pressure,  psf. 

Bow  acceleration,  ft/sq  sec. 

C.G.  acceleration,  ft/sg  sec. 

Fan  power,  horsepower. 

Roll  angle,  degrees. 

Yaw  angle,  degrees. 

Lateral  acceleration  ft/sq  sec. 

Speed  through  the  water,  knots. 

Turn  radius,  feet. 

Plenum  volume,  cu  ft. 

X  displacement,  feet. 

Y  displacement,  feet. 

Air  flow  rate  in,  cb  ft/sec. 

Air  flow  rate  out,,  cb  ft/sec. 

Net  force  in  x  direction,  lb* 

Wave  force  in  x  direction,  lb. 

Starboard  thrust,  lb. 

Port  thrust,  lb. 

Pitch  rate,  deg/sec. 

Roll  rate,  deg/sec. 

Yaw  rate,  deg/sec. 

Rudder  angle,  degrees. 


Number 

Name 

01 

TIME 

02 

ETA 

03 

Z 

04 

THETA 

05 

PB 

06 

BOKACC 

07 

ACC 

08 

FANPWR 

09 

PHI 

10 

BETAS 

11 

ACCLAT 

12 

U 

13 

TRADIUS 

14 

VOLP 

15 

X 

16 

Y 

17 

QIN 

18 

QOUT 

19 

GFXXX 

20 

FXPWAV 

21 

THSTS  (1) 

22 

THSTP  (1) 

23 

QDEG 

24 

PDEG 

25 

RDEG 

26 

DELES 
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LISTING  OF  MAIN  PROGRAM  AND  SUBROUTINES 


Main  Program 

Discussion  -  The  main  program  contains  the  logic  for 
interconnecting  the  various  subroutines.  In  addition,  it 
compares  the  running  value  of  time  tfith  the  finish  time, 
calculates  the  next  value  of  time  for  printing  after  each 
print  time,  etc.  The  main  program  also  contains  the  logic 
for  the  fixed  coordinate  trajectory  (x,  y,  lj)  )  by  the 
trapazoidal  rule. 

Output  -  The  output  from  the  main  program  is  controlled  by 
the  print  option  switch  ITRAJ.  If  the  value  of  this  switch 
is  1,  the  program  will  print  the  current  values  of  the  time 
and  the  craft  velocities  and  displacements  (angular  and 
translational) . 

BLOCK  DATA 

Discussion  -  The  block  data  aie  part  of  the  main  program  and 
contain  all  COHilON-statements  used  in  various  subroutine 
programs.  All  values  therein  are  initialized  with  "0.0" 
before  the  first  data  are  read  into  the  program.  They  keep 
their  current  values  between  consecutive  runs  (controlled  by 
card  0  13)  as  far  as  they  are  not  recalculated  due  to  a 
change  of  input  values  or  within  the  subroutines  themselves. 
Various  constants,  too,  are  transferred  by 
COMMON-statements. 

Subroutine  INCON 

Discussion  -  Subroutine  INCON  contains  the  logic  for  the 
reading  of  all  input  data,  the  initialization   of   variables 
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and  for  the  initiation  of  new  cases.  The  subroutine  is 
programmed  so  that  the  input  data  is  non-structured  i.e. 
there  is  no  set  order  of  input  cards  from  one  read  statement 
to  another.  However,  a  given  read  statement  may  have  more 
than  one  input  card.  In  runs  with  multiple  cases,  only  those 
input  parameters  which  are  to  be  changed  for  the  new  case 
need  be  read  in. 


Usa_ge 


CALL  INCON 


Output  -  The  output  from  INCON  is  discussed  earlier. 

Subroutine  DMINV 

Discussion  -  Subroutine  DI1INV  is  a  matrix  inversion  package 
used  for  inverting  the  mass  matrix.  The  method  used  is  the 
standard  Gauss-Jordan  technigue.  The  determinant  of  the 
matrix  is  also  calculated.  If  the  determinant  is  zero  the 
matrix  is  singular  and  the  program  will  stop. 

Usage  -  CALL  DMINV  (A,  N,  D) 

inscription  cf  Parameters 
A  =  The  input  matrix,  destroyed  in  computation  and  replaced 

by  resultant  inverse 
N  =  The  order  of  matrix  (maximum  =  6) 
D  =  The  resultant  determinant  of  the  matrix 

Subroutine  WAVES 

Discussion  -  Subroutine  HAVES  calculates  the  wave  forces  and 
moments  acting  on  the  craft.  It  also  generates  the  wave 
heights  at  the  various  stations  around  the  seals  and 
sidewalls,  as  well  as  the  bubble  volume  lost  due  to  wave 
elevation.   The  sea  state  for  irregular  seas  is  computed   in 
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WAVES  by  adding  together  a  series  of  regular  waves  with  an 
appropriate  distribution  of  amplitude  and  frequency. 

Usage  -  CALL  WAVES 

Output  -  The  output  from  WAVES  is  controlled  by  print  option 
switch  IWAVES.  When  this  switch  is  set  equal  to  1 ,  tfAVES 
will  print  the  wave  elevation  around  the  seals  and 
sidewalls,  the  wave  elevation  at  the  center  of  gravity  and 
the  volume  reduction  of  the  bubble  plenum  due  to  waves  as 
well  as  the  total  forces  and  moments  acting  on  the  craft  due 
to  waves. 

Subroutine  SiDEWL 

Discussion  -  Subroutine  SIDEWL  calculates  the  forces  and 
moments  acting  on  the  craft  due  to  the  sidewalls  as  well  as 
the  leakage  flow  rates  associated  with  any  gaps  which  open 
under  the  sidewalls  due  either  to  craft  motion  or  waves. 

Usage  -  CALL  SIDEWL 

Output  -  The  output  from  SIDEWL  is  controlled  by  the  print 
option  switch  ISIDWL.  When  this  switch  is  equal  to  1,  the 
program  will  print  the  sidewall  gaps,  immersion  depths  and 
total  forces  and  moments. 

Subroutine  PROP 

Discussion  -  Subroutine  PROP  is  used  to  calculate  forces  and 
moments  on  the  craft  due  to  the  propulsion  system.  The 
subroutine  has  the  capability  of  uneven  thrusting  and 
receives  this  information  by  using  the  function  program  FG1 
to  linearly  interpolate  a  map  input  in  block  16. 
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Usa^e  -  CALL  PROP 

Output  -  The  output  from  PROP  is  controlled  by  print  option 
switch  IPROP.  If  this  switch  is  1,  the  program  will  print 
the  calculated  forces  and  moments. 

Subroutine  RUDDER 

Discussion  -  subroutine  RUDD23  contains  the  logic  for 
calculating  rudder  forces  and  moments  as  well  as  rudder 
motions.  Rudder  motions  are  calculated  by  employing  the 
function  program  FG^  to  linearly  interpolate  an  optional  map 
input. 

Usa^e  -  CALL  RUDDER 

Output  -  The  output  from  rudder  is  controlled  by  print 
option  switch  IRUD.  If  this  switch  is  1,  the  program  will 
print  out  the  total  forces  and  moments  due  to  the  rudder  . 

Subroutine  AEROD 

Discussion  -Subroutine  AEROD  calculates  the  aerodynamic 
forces  and  moments  on  the  craft.  The  subroutine  currently 
uses  curve  fits  of  the  aerodynamic  coefficients  vs  sideslip 
angle. 

Usage  -  Call  AEROD 

Output  -  The  output  from  AEROD  is  controlled  by  print  option 
switch  IAEROD.  If  this  switch  is  1,  the  program  will  print 
the  total  forces  and  moments. 
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Subrourine  INTGRL 

JDi^^Mssion  -  Subroutine  INTGRL  is  used  to  integrate  a  system 
of  first  order  ordinary  differential  equations.  It  uses  a 
variable  time  step  technique  based  on  the  Runge-Kut ta-Merson 
algorithm.  A  fixed  time  step  Runge-Kutta  algorithm  is  also 
available.  The  program  will  stop  the  calculation  if  the 
time  step  becomes  smaller  than  10-6sec. 

Usa^e  -  CALL  INTGRL  (TIME) 

Dsscriotion  cf  Parameters 

TIME  =   current  value  of  the  independent  variable 

Subroutine  RHS 

discussion  -  Subroutine  RHS  is  the  subprogram  containing  the 
FORTRAN  expressions  for  the  right  hand  side  of  the  system  of 
first  order  differential  equations.  It  contains  the  logic  to 
calculate  the  bubble  volume,  area  and  pressure  as  a  function 
of  time  as  well  as  the  logic  which  sums  the  individual 
forces  ,  moments  and  leakages  from  the  various  craft 
components  to  get  craft  totals.  RHS  also  contains  the 
statements  used  for  the  writing  of  the  scratch  file  fcr  the 
tabular  output. 

Usa^e  -  CALL  RHS  (VALUE) 

Description  of  Parameters 

VALUE  =  The  array  containing  the  values  of  the  right  hand 
side  of  the  differential  equations  as  calculated  in 
RHS,  i.e.,  for  a  system  of  equations  of  this  form 
y(i)  =  f(y(i),t)  the  VALUE  array  is  used   to   store 
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the  values  of  f(y(i),t). 

Output  -  The  output  from  RHS  is  controlled  by  print  option 
switch  IRHS.  When  this  switch  is  1,  the  subroutine  will 
print  all  right  hand  side  forces,  moments,.  accelerations, 
pressure,  air  flow  rates  and  fan  power  besides  the  complete 
VAL-  and  VALUE-  arrays.  This  print  switch  is  the  roost 
efficient  one  for  debugging  purposes  but  also  computer  time 
consuming. 

Subroutine  BOtfSL 

Pi^cussion  -  Subroutine  BOWSL  calculates  the  forces  and 
moments  due  to  the  bow  seal  and  the  leakage  flow  rate 
associated  with  any  gaps  which  open  up  under  various  seal 
stations. 

Usa3£  ~  CALL  ROWSL 

Output  -  The  output  from  EOWSL  is  controlled  by  print  option 
switch  IBOWSL.  If  this  switch  is  1,  the  program  will  print 
the  gaps  and  wetted  lengths  at  the  various  stations  on  the 
bowseal  as  well  as  the  total  forces  and  moments  due  to  the 
bowseal. 

Subroutine  STNSL 

Discussion  -  Subroutine   STNSL  calculates   the   forces   and 

moments  acting  on  the  craft  due  to  the  stern  seal  as  well  as 

the  leakage  flow  rates  arising  from  any  gaps  which  open 
under  the  seal. 

Usa^e  -  CALL  STNSL 

Cutput  -  The  output  from  STNSL  is  controlled  by  print  option 
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switch  ISTNSL.  When  this  switch  is  equal  to  1,  the  program 
will  print  the  seal  wetted  length,  and  gaps  as  a  function  of 
station  as  well  as  the  craft  forces  and  moments  due  to  the 
stern  seal. 

Subroutine  FAN 

Discussion  -  Subroutine  FAN  is  used  to  calculate  the  inflow 
to  the  plenum,  stern  seal,  and  bowseal  from  the  fans  for  a 
given  pressure  differential  across  the  fans,  the  subroutine 
uses  tabular  data  and  allows  for  changes  in  fan  speed  as 
well  as  the  number  of  fans. 

Usa^  -  CALL  IAN 

Output  -  The  output  from  FAN  is  controlled  by  print  option 
switch  IRHS.  If  this  switch  is  1 ,  the  program  will  print  the 
flow  rates  and  pressure  differentials  across  the  bowseal, 
cushion  and  stern  seal  fans,  as  well  as  ideal  fan  power, 
actual  fan  power  and  fan  efficiency. 

Subroutine  COLFIL 

Discussion  -  Subroutine  COLFIL  is  used  to  format  the  program 
output.  It  can  provide  up  to  16  variables  listed  in  two 
tabular  summaries  of  8  variables  each.  The  subroutine  has  a 
plotting  package  that  permits  use  of  either  the  print  plot 
or  calcomp  (continuous  line)  option.  Multiple  curves  may  be 
plotted  and  up  to  10  graphs  may  be  drawn  on  each  run. 

Usage  -CALL  COLFIL 

Output  -  The  output  from  COLFIL  is  discussed  earlier. 
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Subroutine  SAM 

Discussion  -  subroutine  SAM  is  supposed  to  calculate  the 
shear  and  bending  moments,  but  is  not  used  in  this  program 
version  and  replaced  by  a  dummy  subroutine. 

Function  FG1 

discussion  -  The  purpose  of  FG1  is  to  perform  the  task  of 
evaluating  a  tabular  function  y  =  f (x) .  the  procedure  used 
is  linear  interpolation  between  data  points.  If  the  input 
value  of  the  independent  variable  does  not  fall  in  the  range 
of  its  input  tabular  values  the  appropriate  table  extreme  is 
held.  The  function  y  must  be  single  valued  and  the  tabular 
values  for  the  independent  variable  must  be  stored  in 
ascending  order  i.e.  XTAB  (J)  <XTAB  (J+1)  .  The  increment 
between  tabular  values  need  nor  be  constant. 

Usage  -Y  =  FG1  (X,N,XTAB,  YTAB,  IMEM) 

P.ggcr i£tion  of  Parameters 
Y     =   The  output  value  of  the  tabular  function 
X     =   The  value  of  the  independent  variable  to  be  used  to 

evaluate  the  function  y  =  f  (x) 
N     =   The  number  of  elements  in  the  XTAB  and  YTAB  arrays 
XTAB  =   The  name  of  the  array  that  contains  the   values   of 

the  independent  variable  in  ascending  order 
YTAB  =   The  name  of  the  array  that  contains  the   values   of 

the   dependent  variable  corresponding  to  the  values 

of  the  independent  variable  stored  in  XTAB 
IMEM  =   Dummy  variable  used  to  store  location  of  last  entry 

into  table,  to  be  used  as  initial  point   of   search 

at  next  entry 
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Function  SHXYAX 

Discussion  -  Function  SHXYAX  calculates  the  shift  of  the 
x-coordinate  of  a  given  location  due  to  rotation  about  tne 
y-axis  in  the  local  reference  frame.  It  is  used  to  attain 
corrected  levers  the  moment  calculations. 

Usage  -  XC=SHXYAX(X,Z,ANGYAX,PI) 

Description  cf  Parameters 

XC      =  Corrected  x-coordinate 

X       =  Original  x-coordinate  of  given  location 

Z       =  Original  z-coordinate  of  given  location 

ANGYAX  =  Angle  of  rotation  about  y-axis 

PI     =  Constant  pi 

Functions  T1  and  T2 

Discussion  -  Functions  T1  and  T2  are  two  small 
FUNCTION-subprograms  used  to  calculate  various  tr ignometric 
relations  used  by  subroutine  WAVES. 
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LISTING  OF  BLOCKS 


The  names  mentioned  in  the  following  summary  are  used  in 
subroutine  INCON.  Not  all  of  them  are  transferred  as  read  to 
other  subprograms.  An  investigation  of  the  corresponding 
COMMON-statement  may  give  further  hints,  but  care  must  be 
taken  as  those  names  may  be  not  identical  in  all 
subroutines . 


Block_1__Executive_  Input 


Option  tag  1  =  Time  steps 
Card  Column  Format   Name 


1    1-5 
6-15 

16-25 
26-35 

36-45 
46-55 


F10    STIKE 


F10 
F10 

F10 

F10 


FTIME 
DELO 

DELPNT 
TPRINO 


Entry 

00101 

Start    time       (generelly      zero)  , 

sec. 

Finish    time,    sec. 

Initial  integration  time  step, 

sec. 

Print  interval,  sec. 

Print  start  time,  sec. 


Option  tag  2  =  Print  switches 

Card  Column  Format   Name       Entry 

1  1-5      I  00102 

2  5        I      IACCEL     Print    switch    for    lateral 

acceleration   and   integration 

time  step. 
10       I      IVEL       Not  used  at  present. 
15       I      ITSAJ      Print  switch  for  main  program. 
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Card  Column  Format   Name      H±I±%±.1    (cont) 

2    20       I      ISIDWL     Print   switch   for   subroutine 

SIDEWL. 
25       I      IBOWSL     Print   switch   for   subroutine 

BOWSL. 
30       I      ISTNSL     Print   switch   for   subroutine 

STNSL. 
35       I      I WAVES     Print   switch   for   subroutine 

WAVES. 
40       I      IRUD       Print   switch   for   subroutine 

RUDDER. 
45       I      IPROP      Print   switch   for   subroutine 

PROP. 
50       I      IAERDD     Print   switch   for   subroutine 

AEROD. 
55       I      IRHS       Print   switch   for   subroutine 

P.HS  and  FAN. 

Option  tag  3  =  Integrator  tolerances 
Card  Column  Format   Name       Entry 

1  1-5      I  00103 

2  1-2      12     NEQS       Number  of   integrators   to   be 

used  (maximum  of  15)  . 
3-4      12     JQQ        Control   tag   for   integration 

procedure. 

00  =  original  Runge-Kutta 
procedure. 

01  =  fixed  step  size 
Runge-Kutta  procedure. 

02  =  variable  step  size 
Runge-Kutta  procedure.  Only 
the  minimum  step  size  is 
output. 
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Card  Column  Format   Name 
3,4,.  1-80     8F10   TOL (J) 


Entry  (cont) 

Integrator   error    tolerance, 
one  per  integrator. 


Note.1 

The  usual  procedure  for  selecting  integrator  tolerances  is 
to  choose  values  on  the  order  of  10-5  to  10-8.  For  a  new 
craft  the  follov.'ing  procedure  is  used: 

Make  a  run  with  a  given  set  of  tolerances;  halve  the 
tolerances  and  rerun.  If  the  results  do  not  cnange 
drastically,  try  doubling  the  original  tolerances  and  rerun. 
As  changes  are  noted  in  the  integrator  outputs  (i.e.  craft 
velocities,  etc.),  hold  the  tolerances  for  those  integrators 
fixed  and  vary  the  others  until  the  largest  tolerances 
possible  are  obtained.  Once  chosen  for  a  given  craft  the 
tolerances  rarely  have  to  be  changed. 


Option  tag  4  =  Tabular  summary  of  output 
Card  Column  Format   Name       Sntry 


1    1-5 


IVERI 


10 

I 

ILATRL 

15 

I 

NVD 

20 

I 

NVI 

25 

I 

NLD 

30 

I 

NLI 

00104 

Print  switch  for  summary  one 
1  =  print,  0  =  no  print. 
Print  switch  for  summary  two 
1  =  print,  0  =  no  print. 
Not  used  at  present. 
Not  used  at  present. 
Not  used  at  present. 
Not  used  at  present. 


Option  tag  5  =  Program  option  switches 
Card  Column  Format   Name       Entry 
1    1-5      I  00105 
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Card  Column  Format   Name       Entry,  (cont) 

1    6-15     F10    I3D0F      Switch    for    lateral    plane 

motions. 
1  =  only  lateral  plane  motions 

0  =  six  degrees  of  freedom. 
16-25    F10    ISRGE      Switch  for  constant  speed. 

1  =  surge  equation  not  used. 

0  =  surge  velocity  allowed  to 
vary. 

26-35    F10    ITRIM      Switch  for  trim. 

1  =  thrust  varied  to  give 
proper  speed.  If  this  option 
is  selected,  then  there  must 
be  input  in  Block  16  an 
estimate  of  thrust  and  a  0.0 
for  the  number  of  data  points. 
0  =  thrust  as  input. 

36-45    F10    IDIA       Switch   for    membrane    study 

concerning  plenum  air  bubble. 

Option  tag  6  =  Program  modification  settings 
Card  Column  Format   Name      JLntry; 
11-5      I  00106 

6-15     F10    PR0M01     Program  modification  setting  1 

(temporarily  used  as  deadrise 
factor) . 

16-25    F10    PPOM02     Program  modification  setting  2 

(temporarily  used  as  roll 
damping  factor) . 

26-35    F10    PHOM03     Program  modification  setting  3 

(temporarily  used  as  print 
switch  for  vertical  forces  and 
roll  moments  in  subroutine 
SIDEKL) . 
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Card  Column  Format   Name 


1    36-45 

F10 

PR0MC4 

46-55 

F10 

PR0M05 

56-65 

F10 

PR0M06 

66-75 

F10 

PRO  11 07 

Entry  (conr) 

Program  modification  setting  4 

Program  modification  setting  5 

Program  modification  setting  6 

Program  modification  setting  7 


Block  2  Mass  and  .inertial  Input 

Option  tag  1  =1   Summary  mass  properties  card 

Not®!    In   this   option   shears   and   moments   will   not  be 

calculated  correctly. 


Card  Column  Format  Name 

1    1-5  I 

6-15  F10  WEIGHT 

16-25  F10  XS 

26-35  F10  ZS 

36-45  F10  AIXX 

46-55  F10  AIYY 

56-65  F10  AIZZ 

66-75  F10  AIXZ 


Entry 

00  201 

Total  craft  weight,  pounds. 

Longitudinal     center     of 

gravity,    feet    forward   of 

transom. 

Vertical   center   of   gravity, 

feet  above  baseline  (keel) . 

Mass  moment  of  inertia  about  x 

axis,  slug-sq  ft. 

Mass  moment  of  inertia  about  y 

axis,  slug-sg  ft. 

Mass  moment  of  inertia  about- z 

axis,  slug-sq  ft. 

Mass   moment   of  inertia  about 

xz  axis,  slug-sq  ft. 


Option  tag  2  =  Discrete  mass  distribution  cards 
Card  Column  Format   Name       Entrv 


1    1-5 


00202 


2,3,.  1-10     F10    AMI  (I)     Discrete  weight,  pounds 
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Card  Column  Format   Name       Entry  (cont) 

2,3,.  11-20    F10    XI(I)      Longitudinal      center     of 

gravity,    feet    forward   of 
transom. 
21-30    F10    YI (I)      Transverse   center  of  gravity, 

feet  to  starboard. 
31-40    F10    ZI  (I)      Vertical   center  of   gravity, 

feet  above  baseline  (keel). 
Last   1-10     F10  Discrete   weights   termination 

card  tag. 

(must  be  less  than  0.0) 

Note^ 

The  maximum  number  of  discrete  weights  is  201.  They  are 
assumed  in  symmetric  distribution  about  the  longitudinal 
center  plane,  so  that  only  weights  on  starboard  side  should 
be  specified,  with  atransverse  location  greater  than  zero. 

Bio ck  3  Craft _Geo m e t r ic ,_I n d u t 

No  options 

Card  Column  Format   Name       Entry 
1    1-5      I  00301 

6-15     F10    NSTA(1)    Number   of   stations   on   port 

sidewall. 
16-25    F10    NSTA(2)    Number     of     stations    on 

starboard  sidewall. 
26-35    F10    NSTA(3)    Number   of   stations   on    bow 

seal. 
36-45    F10    NSTA(4)    Number  of   stations   on   stern 

seal. 
46-55    F10    XLTOT      Total  craft  length. 
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Kotej. 

The  maximum  number  of  stations  is  11  for  each  element.  The 
program  will  automatically  subdivide  the  element  into  N- 1 
sections  (N  is  the  number  of  stations)  and  calculate  the  x 
and  y  coordinates  of  each  station. 

Block  4  Sidewall  Input 

Option  tag  1  =  Sidewall  input 
£§£^  Column  Format   Nam?       Entry 
1    1-5      I  00401 

6-15    F10    YSV;       y   distance  from  centerline  to 

sidewall,  feet. 
16-25    F10    XLSH       Average    wetted    length    of 

sidewall,  feet. 
26-35    F10    CFSK      Leakage   orfice  coefficient  of 

sidewall. 
36-45    F10    CDSW       Cross-flow  drag  coefficient  of 

sidewall . 
46-55    F10    &VBASW     Average    beam    of   sidewall, 

feet. 

Notej. 

Under  this  option  tag  the  program  will  also  read  the 
sidewall  integral  table  file  stored  on  a  special  unit  for 
fast  access. 

2l2SJS_5_Stern_Seal_In2Ut 

No  options 

£i*£d  Column  Format   Name       Entry, 

1    1-5      I  00501 


88 


Card  Column  Format   Name       Entry  (cont) 

1    6-15     F10    XSSI       x  coordinate   of   seal   hinge, 

feet  forward  of  transom. 
16-25    F10    ZSSI       z  coordinate   of   seal   hinge, 

feet  above  keel. 
26-35    F10    ALEAK      Base  leakage  area,  sg  ft. 

(3.79  sg  ft  is  nominal.) 
36-45    F10    CFSS      Seal      leakage      orifice 

coefficient . 
46-55    F10    ELMAXS     Length  of  rear  support   cable, 

feet. 

(1.875  feet  is  full  down.) 
56-65    F10    DPBS       Pressure   differential  between 

stern  seal  and  plenum  chamber, 

psf . 
66-75    F10    XLF        Length  of  leading  edge  of   the 

seal,  feet. 

Block_6 ..Bcwseal  Input 

Option  tag  1  =  First  data  card 
Card  Column.  Format   Name       Entry 
1    1-5      I  00601 

6-15     F10    XBSI       x  coordinate   of   seal   hinge, 

feet  foreard  of  transom. 
16-25    F10    CFBS       Seal       leakage       orifice 

coefficient . 
26-35    F10    DPBS       Pressure  differential   between 

bow   seal   and  plenum  chamber, 

psf . 
36-45    F10    ZBSI       z   coordinate   of   seal  hinge, 

feet  above  keel. 
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Card   Column    Format      Nam®  ILlii£X    (cont) 

1        46-55        F10         ELMAXB  Length      of   rear   support    cable, 

feet. 

(1.875  feet  is  full  down.) 
56-65    F10    XBF        Length  of  seal   leading   edge, 

feet. 
66-75    F10    BLEAK     Base  leakage  area,  sg  ft. 

Option  tag  2  =  Second  data  card 

Card  Column  Format   Name      Entry 

1    1-5      I  00602 

6-15     F10    CEKCAB     Length    of    middle    support 

cable,  feet. 
(1.875  feet  is  full  down.) 

Block  7  Plenum  Input 

Option  tag  1  =  Plenum  geometry 
Card  Column  Format   Name       HQ££I 
I  00701 

6-15     F10    XLBW       Plenum    length     at     water 

surface,  feet. 
16-25    F10    X5BW       Plenum  width  at  water  surface, 

feet. 
26-35    F10    XPWV       x  coordinate  of  pressure   wave 

pivot   point,   feet  forward  of 

transom. 
36-45    F10    WIDTH      Plenum  width  at  deck,  feet. 
46-55    F10    XL        Plenum  length  at  deck,  feet. 
56-65    F10    XCPO       x    coordinate   of   center   of 

pressure,   feet    forward   of 

transom. 
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Card   Column    r oroiat      Name  Entry_    (cont) 

1         66-75         F10         BUBHGT  Plenum    average    height,    feet. 

Option    tag    2    =   Critical    Froude    number 
Card    Column    Format      Name  Entry 


1 


1-5 


00702 


6-15     F10    FNCRIT     Froude  number  corresponding  to 

hump  speed. 

ll2£k_3_Propulsion_Input 


No  options 

Card  Column  Format   Name 

1    1-5  I 

6-15  F10    XPO 


16-25    F10    YPO 


26-35    F10    ZPO 


Entry 

00801 

x    coordinate    of   propeller 

center,    feet    forward    of 

transom. 

y  distance  from  centerline   to 

propeller  center,  feet. 

z   coordinate   of    propeller, 

feet  above  keel. 


Block  9  Rudder  Incut 


No  options 

Card  Column  Format   Name 


1 


1-5 


6-15 


F10 


XRO 


Entry 


00901 


x  coordinate  of  centroid  of 
rudder,  feet  forward  of 
transom. 
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Card  Column  For  mat   Name 
1    16-25    F10    Yfi 

26-35    F10    ZRO 


36- 

-45 

F10 

RSPAK 

46- 

-55 

F10 

RASPR 

56- 

-65 

F10 

RAF.EA 

66- 

-75 

F10 

RTC 

Entry  (cont) 

y   distance  from  centerline  to 

rudder  centroid,  feet. 

z   coordinate   of   centroid  of 

rudder,  feet  above  keel. 

Rudder  span,  feet. 

Rudder  aspect  ratio. 

Rudder  area,  sq  ft. 

Average   thickness   ratio   of 

rudder  section. 


Block  _10_Aerodvj^amic_In£ut 

No  options 

Card  Column  Format   Name 


1 


1-5 


Entry 
01001 


6-15     F10    XLAERO  Reference  length,  feet. 

16-25    F10    BEAM  Reference  width,  feet. 

Notex 

The   reference   length   and  width    are    those    used   to 

nondimmensionalize  the  wind  tunnel  data  used  in  the  program. 

Block, 11  Sea  State  Input 


Option  tag  1  =  Amplitude  versus  frequency 
Card  Column  Format   Name       Entry 
1    1-5      I  01101 


6-15 


F10 


16-25    F10 


NWAVE 


BETAD 


Number   of   wave   components, 
maximum  of  10. 

Initial  wave  heading,  degrees. 
(180  degrees  =  ahead  waves) 
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Card  Column  Format   Name 


Entry;  (cont) 


2,3,.  1-10     F10    OMEGA (I) 
11-20    F10    AW (I) 


Wave  frequency,  rad/sec, 
Wave  amplitude,  feet. 


Option  tag  2  =  Amplitude  versus  wave  length 
Card  Column  Format   Name      HEt£Z 
1    1-5      I  01102 


6-15 


F10    NHAVE 


16-25         F10         BETAD 


2,3,.  1-10     F10 
11-20    F10 


Number   of   wave    components, 

maximum  of  10. 

Initial  wave  heading,  degrees. 

(180  degrees  =  ahead  waves) 
WAVLEN(I)  Wave  length,  feet. 
Atat(I)      Wave  amplitude,  feet. 


Note: 

In  this  options  irregular  seas  are  constructed  by  addition 
of  regular  wave  components.  For  a  regular  sea  case  only  one 
component  is  input.  If  the  number  of  components  is  zero 
(calm  water) ,  only  the  first  card  may  appear. 

Option  tag  3  =  Wave  components  by  calculation  for  given 

period  limits 


Card  Column  Format  Name 
1    1-5      I 

6-15     F10  NWAVE 

16-25    F10  BETAD 

26-35    F10  SHTWV 

36-45    F10  PERL 


Entry 

01103 

Number    of   wave   components, 

maximum  of  10. 

Initial  wave  heading,  degrees. 

(180  degrees  =  ahead  waves) 

Average   height   of1/3  highest 

waves,  feet. 

Shortest  significant  period  of 

waves,  seconds. 
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Card  Colaran  Format   Name       Ent_ry  (cont) 

1    46-55    F10    PERH       Longest   significant  period  pf 

waves,  seconds. 

Option  tag  4  =  Wave  components  by  calculation  for  given 

frequency  limits 

Card  Column  Format   Name       Entry. 

1    1-5      I  01104 

6-15     F10    NWAVE      Number    of   wave   components, 

maximum  of  10. 
16-25    F10    BETAD      Initial  wave  heading,  degrees. 

(180  degrees  =  ahead  waves) 
26-35    F10    SHTWV      Average   height   of1/3  highest 

waves,  feet. 
36-45    F10    WWN        Lowest   significant   frequency 

of  waves,  Hertz. 
46-55    F10    WWX        Highest   significant  frequency 

of  waves,  Hertz. 

Notej. 

In  options  3  and  4  the  wave  components  are  calculated  in 
order  to  ease  the  input  requirements  for  irregular  sea 
cases. 

Block  12  Initial  Conditions 

No  options 

Card  Column  Format   Name       Entry. 

1    1-5      I  01201 

6-15     F10    UO        Initial  speed,  knots. 
16-25    F10    THETO      Initial   pitch  angle  (positive 

-  nose  up,  in  degrees) . 
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Card  Column  Format   Name 


1    26-35    F10    DSO 


36-45    F10    DELPI 


Entry,  (cont) 

Initial  draft   at   center   of 

gravity,  inches. 

Initial  bubble  gage  pressure, 
psf . 


Block_13_End_of__Case 

No  options 

Card    Column    Format      Name  Entry. 

1         1-5  I  013 

Notej^ 

This   card   signals   the   end   of  data  for  a  given  case.  Two 

consecutive  013  cards  with  no  data  between  will  stop  the  job 

unduly. 

Block    14    End   of   Run 

No    options 

Card  Column  Format   Name       Entry 

1    1-5      I  014 

Note^_ 

This  card  signals  the  end  of  the  run.  this  will  be  indicated 

by  the  message:  COMPLETED  ALL  RUNS. 

Block _16  Thrust  and  Rudder  Input 

Option  tag  1  =  Starboard  thrust  input 
C.ard  Column  Format   Name       Entry 
1    1-5      I  01601 
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Card  Column  Format   Name       Entrv_  (cont) 

1    6-15     F10    THST1      Value  of  the   thrust   for   the 

starboard  screw.  This  is  an 
estimate  or  kept  fixed 
depending  on  option  switch  in 
Block  00105. 

16-25    F10    NPS        Number   of   data  points  in  the 

thrust  nap.  Left  blank  if  no 
thrust  map  is  used. 

26-35    F10    STHS       Side  thrust  coefficient. 
2,3,.  1-80     8F10   TIS (J)     Data   points   for   independent 

variable  time  in  ascending 
order. 

1-80     8F10   THSTS(J)   Corresponding  data  points   for 

the  dependent  variable. 

lots!   NPS  must  not  be  larger  than  25. 

Option  tag  2  =  Port  thrust  input 

Card  Column  Format   Name       Entry 

1    1-5      I  01602 

6-15  F10  THST2  As  above  for  port  screw. 

16-25  F10  NPP  As  above  for  port  screw. 

26-35  F10  STHP  As  above  for  port  screw. 

2,3,.  1-80  8F10  TIP  (J)  As  above  for  port  screw. 

1-80  8F10  THSTP(J)  As  above  for  port  screw. 

No^fLL   NPP  must  not  be  larger  than  25. 

Option  tag  3  =  Rudder  motion  input 
Card  Column  Format   Name       Entry 
1    1-5      I  01603 
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Card  Column  Format   Name 
1    6-15     F10    DELP 


Entry_  (cont) 


16-25    F10    NPR 


2,3,.  1-80     8F10   TIP  (J) 


1-80 


Value  of  rudder  if  constant 
rudder  is  desired  (in 
degrees) . 

Number  of  data  points   of   the 
rudder    map.   Left   blank   if 
constant  rudder  is  desired. 
Values  of  independent  variable 
time  in  ascending  order. 
F10  DELRUD (J)   Corresponding  values  of  rudder 

angle  (in  degrees) . 
No_t.e_j^    Positive   rudder  angle 
is  rignt  rudder. 


Note:   NPR  must  not  be  larger  than  25. 


Block  17  Bow  and  Stern  Seal  Pressure  Differences  Input 


Option  tag  1  =  Bow  seal 
Card  Column  Format   Name 
1    1-5      J 

6-15     F10    NB 


2,3,.  1-30 
6,7,.  1-80 


8F10   TKEB(I) 
8F10   DELB(I) 


Entry 

01701 

Number   of  data  points  for  the 

bow  seal. 

Time    points    in    ascending 

order. 

Corresponding   differences   in 

bow  seal  pressures. 


Option  tag  2  =  Stern  seal 
Card  Column  Format   Name       Entry, 
1    1-5      I  01702 
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Card  Column  Format   Name 


1    6-15     F10    NS 


2,3,.  1-80     8F10   TMES(I) 


6,7,.  1-80     8F10   DETS(I) 


Entry_  (cont) 

Number   of  data  points  for  the 

stern  seal. 

As   for   option   1,   for  stern 

seal. 

As  for  option  1,  for  the  stern 

seal. 


Nptei 

If  constant  pressures  are  desired,  those  values  shall  be 
■input  in  block  5  and/or  block  6.  This  block  can  be  used  for 
one  seal  alone,  too. 

Block  1 8  Title  Card 


No  options 

Card  Column  Format   Name 


1    1-5 


1-80 


20A4   TITLC 


Entr^ 

01801 

Alphanumeric  data  to  be 
printed  as  title  in  input 
summary. 


Block_J_9_Fan__Maps 


Option  tag  1  =  Bow  seal  fans 
Card  Column  Format   Name       Entry_ 
1    1-5      I  01901 


6-15 

F10 

ENBFAN 

16-25 

F10 

BRPM 

26-35 

F10  . 

NPTSB 

Number  of  bow  seal  fans. 
Speed  of  bow  seal  fans,  rpm. 
Number  of  data  points, 
(maximum  of  25) 
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Card  Column  Format   Name 


Entry  (cont) 


1    36-45    F10    READIN 


2,3,.     1-80     8F10   PBFAN(J) 


6,1,.     1-80     8F10   QBFAN(J) 


If  this  is  a  1.0  then  the 
following  data  points  are 
read.  If  left  blank  no  further 
read  takes  place. 
Tabular  values  of  fan  pressure 
differential,  psf. 
Corresponding  volumetric  flow 
rate,  cb  ft/sec. 


Option  tag  2  =  Cushion  fans 
Card  Column  Format   Name       Sntry_ 
1    1-5      I  01902 


6-15 

F10 

ENHFAN 

16-25 

F10 

EMRPK 

26-35 

F10 

NPTSi"! 

36-45 

F10 

READIN 

2,3,.  1-80 
6,7,.  1-80 


8F10   PKFAN(J) 


8F10   QMFAN(J) 


Number  of  cushion  fans. 
Speed  of  cushion  fans,  rpm. 
Number  of  data  points, 
(maximum  of  25) 

If  this  is  a  1.0  then  the 
following  data  points  are 
read.  If  left  blank  no  further 
read  takes  place. 
Tabular  values  of  fan  pressure 
differential,  psf. 
Corresponding  volumetric  flow 
rate,  cb    ft/sec. 


Option  tag  3  =  Stern  seal  fans 

Card  Column  Format   Name       Entry 

1    1-5      I  01903 

6-15     F10    ENSFAN     Number  of  stern  seal  fans. 
16-25    F10    SRPM       Speed  of  stern  seal  fans,  rpm 
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Card  Column  Format   Name 
1    26-35    F10    NPTSS 

36-45    F10    RSADIN 


2,3,.  1-8G     8F10   PSFAN(J) 


6,7,.  1-80     8F10   QSFAN(J) 


Fntr_y_  (cont) 

Number  of  data  points. 

(maximum  of  25) 

If   this   is   a   1.0   then  the 

following   data    points    are 

read.  If  left  blank  no  further 

read  takes  place. 

Tabular  values  of  fan  pressure 

differential,  psf. 

Corresponding   volumetric  flow 

rate,  cb  ft/sec. 


Note^ 

The  tabular  pressure  data  for  all  the  fans  must  be  stored  in 

order  of  increasing  pressure.  The  data  is  for  a  single  fan. 

Block  2  0  Required  Curves __or.  ..Graph 


No  options 

Card  Column  Format   Name 


1    1-3 


1-10 


1011   NCURV 


Entry 

020 

The  number  of  this  curve  on 
the  corresponding  graph,  i.e. 
column  1  pertains  to  graph  1, 
column  2  to  graph  2,  e.t.c. 

0  =  only  one  curve  per  graph. 

1  =  first  of  more  than  one 
curves  per  graph. 

2  =  intermediate  of  three 
curves  per  graph. 

3  =  last  of  more  than  one 
curves  per  graph. 
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NotejL  Maximum  number  of  curves  per  graph  =  3. 
Maximum  number  of  graphs  =  10. 

Block  21  Tabular  Summaries 


No  options 

Card  Column  Format   Name 

1  1-3  13 

2  1-16  812    I5UM1 


1-16     812    ISUM2 


Entry 

021 

Numbers  corresponding  to  the 
variable  names  to  be  printed 
out  in  summary  one. 
Number  corresponding  to  the 
variable  names  to  be  printed 
out  in  summary  two.  Sumraary 
two  is  printed  out  in 
successive  runs  without  the 
cards  being  repeated. 


ElQgJi_22  Re gu ire d__ Graphs 


No  options 

Card    Column    Format      Name 


1    1-3 


1-2 


3-4 


13 


12 


12 


NGRAF 


NDRW 


1-52     2612   NXYS 


Entry 

022 

Number 


of    graphs   desired, 


to   and 


up 


right   justified, 
including  10. 

00  =  print  plot 

01  =  CALCOMP  plot 

Numbers  corresponding  to  the 
variable  names  desired  for  the 
X  and  Y  axes  of  the  graphs  in 
alternating  manner. 
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Card   Column    Format      Name 
4         1-48  6A8         TICRD 


Entrx  (cont) 

Alpha-numerical  title  on  the 
first  line  of  the  CALCOMP 
plot.  Blank  if  print  plot  is 
desired. 


Final_Fote 


in  order  to  maintain  the  pitch  angle  curve  of   figure 
the   present   computer   model   has   a   speed  dependent  shift 
routine  of  XCPO  of  the  following  form: 

XCPU=XCP+A*(U*0. 5921-30.0)  **2+B 
After  program  modifications  the  coefficients  A  and  B  may  be 
re-calculated  by  setting  XCPO  equal  to  XCP  and  using 
different  input  values  for  XCPO  at  10  knots  and  30  Knots 
until  the  desired  pitch  angles  are  obtained.  Figure  8  gxves 
some  initial  estimates.  From  the  two  dara  points  the  valu 
for  A  and  B  can  be  determined  replacing  XCPO  by  its  original 
value  again. 
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PLOT   1 


CC5 


Old 


020 


02E 


030 


■SCRLE=5.00E+00   UNITS    INCH. 
■SCRLE=5.00E-01   UNITS    INCH. 

RQGRnM   2,    1..    STRniGHT   RUN   RT   20   KNOTS 
LOT    IS   PITCH   RNGLE  UERSUS   TIME 
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PLOT   ? 


OiO 


02C 


025 


030 


crRL[- 


juriLc.==5.00E+00  UNITS  INCH 
SCflLE=1.00E+00  UNITS  INCH 
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3QGRRM   2 
-0T    IS   I 


1 
_L 


STRRIGHT   RUN   RT   20  KNO 
DISPLRCEMENT        UERSUS   TIME 


S 


jOLf 


PLOT    3 


005 


OiQ 


01! 


023 


Q2t 


030 


03-' 


T5    INCH. 


-SCRLE=5.00E+Q0  UNI 
-SCALE=5.00E+00  UNITS    INCH. 

TORRM   2,    1,    STRAIGHT    RUN   RT    20   KNOTS 
'LOT    IS   PLENUM   PRESSURE      UERSUS   TIME 
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PLOT   4 


010 


02: 


025 


030 


03T 


SCRLE=5.00E+00  UNITS  INCH= 
SCRLE=5.00E+01  UNITS  INCH. 

RGGRRM  2,  1.  STRRIGHT  RUN  RT  20  KNOTS 
QT  IS  THRUST  STRRBORRD  UERSUS  TINE 
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PLOT  5 


010 


CIS 


020 


025 


n.in 


03r 


-SCRLE=5.00E+00  UNITS  INCH- 
SCRLE=2.00E-01  UNITS  INCH. 

RQGRRM   2,    2=    STRniGHT 
LOT    IS   PITCH   RNGLE 
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RUN   RT    20   KNOTS 
UERSUS   TIME 


PLOT  6 


010  315  020  025  030 


SCALE-5.00E+00  UNITS    INCH.  io8 

SCHLE=1.00E+00   UNITS    INCH. 

RQGRRM   2,    2.    STRRIGHT  RUN   RT   20   KNOTS 
-QT    IS  Z   DISPLACEMENT        UERSUS   TIME 


/pf 


PLOT  7 


s. 


nor 


010 


015 


020 


Vcb 


030 


-SCRLE--5.00E+00   UNITS    INCH. 
-SCfiLF..---5.00E+00  UNITS    INCH. 

TORRM   2,    2*    STRfllGHT   RUN   RT   20   KNOTS 
LOT    IS   PLENUM   PRESSURE      UERSUS   TIME 
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PLOT  8 


010  015  020  "'  025  030 


SCflLE=5.00E+00  UNITS  INCH. 
SCflLE=5.00E+0i  UNITS  INCH. 

RQGRflM  2,  2*  STRAIGHT  RUN  RT  20  KNOTS 
LOT  IS  THRUST  STRRBORRD  UERSUS  TIME 
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PLOT   9 


010 


0i5 


020 


:2b 


030 


035 


1  » 


-3CRLE=5.00E+00  UNITS  INCH, 
-SCRLE=5.00E-01  UNITS  INCH 

ROGRRN 


IS 


3*    1.    STRRIGH 
PITCH   RNC 


bl_t 


RUN   RT   20   KNOTS 
UERSUS   TINE 
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PLOT   10 


OiO 


02C 


02b 


030 


■SCflLE=5.00E+00  UNITS   INCH,  112 

-SCRLE=-\L  OQE+00   UNITS    INCH, 

ROGRRM  3^  1.  STRRIGHT  RUN  RT  20  KNOTS 

LOT  IS  Z   DISPLRCEMENT   UERSUS  TIME 

at 


PLOT  11 


005 


010 


02C 


02b 


030 


SGRLE--5.00E+00  UNITS   INCH. 
-SCRLE-5. 00E+00  UNITS   INCH. 

ROGRRM   3,    1;    STRniGHT   RUN   RT    20   KNOTS 


LOT    IS   PLENUM   PRESSURE      UERSUS   TIME 
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PLOT  12 


OiO  CIS    "  "    020  02b  030  035 


SCRLE=5.QQE+00   UNIT3    INCH. 
■SCRLE-5.Q0E+01   UN1T5    INCH. 

30GRRM   3,    1.    STRRIGHT   RUN   RT   20   KNOTS 
-QT    IS   THRUST   STRRBORRD   UERSUS   TIME 
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PLOT  13 


010 


02C 


025 


030 


SCRLE=5.00E+00  UNITS  INCH, 
SCRLE=1.00E-Q1  UNITS  INCH. 

ROGRRN  3,  2-  STRniGHT  RUN  RT  20  KNOTS 
LOT  IS  PITCH  ANGLE      UERSUS  TINE 


115 


PLOT  14 


010   ~        G15  Q20  025  030  OjI 


5CRLE-5.00E+00  UNITS    INCH. 
3CRLE---1.00E+00   UNITS    INCH. 

10GRRM   3.    2-   STRRIGHT   RUN   RT   20  KNOTS 
-QT    IS   Z   DISPLACEMENT        UERSUS   TIME 


116 


PLOT   15 


OiO  C15  020  C2b  030 


5CRLE=5.00E+00  UNITS    INCH. 
SCRLE-5.0QE+00   UNITS    INCH, 

iOGRRM   3^    2»    STRRIGHT   RUN   RT   20   KNOTS 
-QT    IS   PLENUM   PRESSURE      UERSUS   TIME 


117 


PLOT   16 


010  0:5  020  025  030 


■SCALED.  00E+00  UNITS  INCH. ' 
-SCRLE=--5,00E+01  UNITS  INCH. 

RQGRRM  3,  2-  STRAIGHT  RUN  RT  20  KNOTS 
^OT  IS  THRUST  STRRBORRD  UERSUS  TIME 


118 


PLOT  17 


030 


03" 


-SCfiLE=5.00E+00  UNITS  INCH. 
-SCRLE=5.00E-Q1  UNITS  INCH. 

ROGRRM  2,  TURN  RT  20  KNOTS  WITHOUT  R>D 
LOT  IS  PITCH  RNGLE      UERSUS  TIME 


119 


PLOT  18 


:2b 


030 


■„: 


SCALE=5.00E+00  UNITS  INCH. 
SCflLE=--2.00E+00  UNITS  INCH. 

RQGRRM   2,    TURN   RT   20   KNOTS   WITHOUT  -R.D 
LOT    IS        PITCH   RRTE  UERSUS   TIME 


120 


PLOT   1Q 


OiO 


020 


02E 


030 


03 


SCRLE==5.Q0E+00   UNITS    INCH- 
3CRLE-2.00E+00   UNITS    INCH. 

3QGRFIM   2,    TURN   RT   20   KNOTS   WITHOUT   R»D 
-QT    IS   ROLL   RNGLE  UERSUS   TINE 


121 


PLOT   20 


^Z^ 


020 


02L 


OiO 


SCRLE=5.00E+00  UNITS    INCH 
SCflLE=5.00E+00  UNITS    INCH 


hOGRRM  2, 
-QT  15 


TURN  RT  20  KNOTS  HITHOUT  R»D 
ROLL  RRTE     UERSUS  TINE 


122 


PLOT   21 


OCL 


010 


015 


020 


02Jj 


030 


SCRLE=5.00E+00  UNITS   INCH. 
SCflLE=5.00E-01   UNITS    INCH. 

hOGRRM  3,  TURN  RT  20  KNOTS  WITHOUT  R*D 


-01  IS  PITCH  ANGLE 


UERSUS  TIME 


123 


PLOT   22 


■SCRLE=5.00E+00  UNITS   INCH, 
■SGRLE=1.00E+00  UNITS    INCH. 

3GGRRM  3,  TURN  RT  20  KNOTS  WITHOUT  R*D 
-OT  IS   PITCH  RRTE     UERSUS  TIME 


124 


PLOT   23 


010   ~  0:5  020  025  030 


)CRLE=5.00E+00  UNITS  INCH, 
3dnLE----2-OOE+00  UNITS  INCH, 

IOGRRM  3,  TURN  RT  20  KNOTS . WITHOUT  R-D* 

QT  IS  ROLL  RNGLE       UERSUS  TIME 


125 


PLOT   24 


•u  cXj 


030 


5CRLE-5.00E+00  UNITS    INCH. 
3CRLE-5.00E+00   UNITS    INCH. 

<0GRRM   3,    TURN   RT   20   KNOTS   WITHOUT   R.D 
QT    IS  ROLL   RRTE  UERSUS   TIME 


126 


PLOT   25 


RLE>5.00E+00   UNITS    INCH 
3CRLE---5.00f.-01   UNITS    INCH 


OiO 


02: 


:2r. 


030 


-1QGRRM   2,    TURN   RT   20   KNOTS   WITH   R.  D. 
QT    IS   PITCH   RNGLE  UERSUS   TIME 


1 


127 


PLOT  ?6 


CRLE---5.00E+00  UNITS  INCH, 
CFLF.---1.00F.+00  UNITS  INCH. 

QGRflM  2,    TURN  RT  20  KNOTS  WITH  R.D 
OT  15   PITCH  RATE     UERSU5  TIME 


1 


120 


PLOT   ?7 


OiO 


02C 


One. 


OiO 


3CRLE-5.0CE+00  UNITS   INCH. 
3CALE-=l„00E+00  UNITS    INCH. 

iOGRRM   2,    TURN   RT    20   KNOTS   WITH   R,D 
-0T    IS   ROLL   RNGLE  UERSUS   TIME 


1 


129 


PLOT   28 


3CRLF.-5.00E+00   UNITS    INCH. 
3CRLE-2.00E+Q0  UNITS   INCH. 

i'QGRRM  2,    TURN  RT  20  KNOTS  NITH  R*D 

-OT  IS    ROlL  RRTE    uersus  time 


1 

.L 


130 


PLOT    ?9 


00^ 


010 


015 


02: 


025 


030 


>CRLE=5.00E+Q0   UNITS    INCH. 
1CRLE--2.00E-01   UNITS    INCH. 

;QGRRM  3,  TURN  RT  20  KNOTS  WITH  R.^D 


OT  IS  PITCH  RNGLE 


ERSUS 


TIME 


131 


PLOT  }0 


025 


030 


3CnLE=5.00E+00  UNITS  INCH. 
3CRLF.--i.00t +00  UNITS  INCH. 

}QGRRM  3^  TURN  RT  20  KNOTS  WITH  R.D 
-QT  IS 


1 


PITCH  RRTE 


UERSUS  TIME 


132 


PLOT    31 


OGS 


010 


015 


020 


ca> 


030 


3CRLE=5.00E+00  UNITS  INCH. 
3CRLE-1.00E+00  UNITS  INCH. 

10GRRM   3*    TURN   RT   20   KNOTS   WITH   R»DC 


1 


OT    IS   ROLL   RNGLE 


UERSUS   TIM 


E 


133 


PLOT  32 


02E 


030 


SCRLE-5.00E+00  UNITS  INCH. 
SCRLE=5.00E+00  UNITS  INCH, 

30GRRM  3^  TURN  RT  20  KNOTS  WITH  R>D 
-QT  IS    ROLL  RRTE     UERSUS  TIME 


L 


134 


PLOT   33 


005 


010 


02C 


025 


030 


0^ 


dCRLE=5.00E+00   UNITS    INCH. 
3CRLE=2.QQE-Q1   UNITS    INCH. 

^QGRRM  3,  TURN  RT  20  KNOTS  WITH  R»D» 
QT  IS  PITCH  RNGLE      UERSUS  TIME 


2 


135 


PLOT  34 


SCALERS.  OOE'+OO  UNITS  INCH. 
SCnLE----1.00E+00  UNITS  INCH. 

^OGRRM  3,  TURN  RT  20 
-0T  IS   PITCH  RRTE 


OTS  WITH  R»D»  2 
UERSUS  TIME 


136 


PLOT  35 


010 


Ci& 


02C 


)2o 


030 


>CRLE=5.00E+00  UNITS   INCH. 
)CRLE----l>00E+00  UNITS   INCH. 

1QGRRM   3^    TURN   RT   20   KNOTS   NITH   R*D* 


<L 


OT    IS   ROLL   RNGLE 


UERSUS   TIME 


137 


PLOT  36 


SCRLE=5.00E+00  UNITS  INCH, 
SCRLE-5. 00E+00  UNITS  INCH, 

WGRRM  3-  TURN  AT  20  KNOTS  WITH  R*D 
.OT  IS    ROLL  ROTE     UER3US  flME 


138 


PLOT   37 


OiO 


Ci5 


022 


121 


030 


3CRLE----5.00E+00   UNITS    INCH. 
3CRLE--2.00E-Q1    UNITS    INCH. 

iOGRRM   3*    TURN   RT    20   KNOTS   WITH   R.^D 
-QT    IS   PITCH   RNGLE  UERSUS   TIME 


139 


PLOT   38 


025 


030 


:rle=--5.ooe+oo  units  inch. 
;rle==i.ooe+oo  units  inch. 

J'GRRN   3,    TURN   RT    20   KNOTS   WITH   R.D. 

3T    15        PITCH   RRTE  UERSUS   TINE 


, — v 

J 


1A0 


PLOT  39 


0C5 


010 


CIS 


020 


025 


030 


CRLD--5.00E+00   UNITS    INCH- 
CRLE>1.00E>00   UNITS    INCH. 

QGRRM  3..  TURN  RT 
QT  IS  ROLL  RNGLE 


20  KNOTS  WITH  R.D-  3 
UERSUS  TINE 


141 


PLOT  40 


^c^ 


'020  025  03Q 


SCRLE--=5.00E+00  UNITS  INCH, 
SCRLE--=5.00E+00  UNITS  INCH, 

POGROM   3*    TURN   RT   20   KNOTS   WITH   R»D*    3 
-0T    IS  ROLL   RATE  UERSUS   TIME 


142 


PLOT   41 


I        ' 


I    f 


»    l' 


OCE 


OJO 


020 


030 


or 


»CflLF^5.00E+00   UNITS    INCH. 
iCfiL£--:2.00E-Ql   UNITS    INCH. 

OGRRM  3,  TURN  RT  20  KNOTS  WITH  R.[ 
QT  IS  PITCH  RNGlE      UERSUS  TINE 


I     B 


4 


143 


PLOT  42 


o;r. 


030 


CALE>5.00E+00  UNITS    INCH. 
CHLE--4.00E+00   UNITS    INCH. 

OGRRM   3.    TURN   RT   20   KNOTS   WITH   R=D 
OT    IS        PITCH   RRTE  UERSUS   TINE 


d 


144 


PLOT   43 


o:z 


OiO 


g:5 


in: 


030 


CRLf>--5.00E+00  UNITS    INCH. 

chle^i.ooe:+oo  units  inch. 

QGRRM   3.    TURN   RT   20   KNOTS   WITH   R.D. 

OT    IS   ROuL   RNGLE  UERSUS   TINE 


d 


145 


PLOT  44 


050 


3CflLE>5.00E:+00  UNITS   INCH. 
3CRLE>2.00E>00   UNITS   INCH 


10GRRM 

QT  IS 


Q 


TURN  RT  20  KNOTS' WITH  R.D. 


ROLL  RRTE 


UERSUS  TIME 


146 


PLOT   45 


010 


C2C 


325 


030 


3CRLE>5.00E>00  UNITS    INCH. 
3CRLi>2.00E-01   UNITS    INCH, 


iQGRRM  3,  TURN  RT  20  KNOTS  WITH  R*D 
-0T  15  PITCH  RNGLE      UFRSUS  TINE 


c~ 


147 


PLOT  46 


02: 


o;o 


5CRLE----5.0QE+00  UNITS  INCH, 
3CRLF.-1.00E+00  UNITS  INCH. 

iOGRRM   3<    TURN   RT    20   KNOTS   WITH   R»D 
-OT    IS        PITCH   RRTE  UERSUS   TIME 


148 


PLOT  47 


V\/> 


OiO 


02^. 


C?L 


OoQ 


3CRLE-5.00[>00  UNITS   INCH, 
3CRLE-.1..00E+00  UNITS    INCH, 

}0GRRM   3*    TURN   RT 


G 


r 


is 


ROLL   RNGLE 


20   KNOTS   NITH   R.D 
ULRSUS   TIME 


1-19 


PLOT   48 


4/VWw^ 


"C2J  C2G  OoO 


3CRLE-5.00E+00  UNITS   INCH- 
3CRLF.==5- 00E+00  UNITS   INCH. 

iOGRRM   3-    TURN   RT   20   KNOTS   WITH   R.D-    E 
.OT    15  ROLL   RRTE  UERSUS   TIME 


150 


PLOT  49 


0-0 


02C 


:?o 


OoO 


scrle^.qoe+oo  units  inch, 
sc;rle>--2.oof:-  01  units  inch. 


3QGRFIM 


3 


TURN   Rl    20 


OT  IS  PITCH  ANGLE 


WOTS  WITH  R.D.  6 
UERSUS  TIME 


151 


PLOT  50 


J2G 


3CflLF.----5.00E:+00  UNITS  INCH. 
iCRLE:---1.00E+00  UNITS  INCH. 

OGRflM  3,  TURN  RT  20  KNOTS  WITH  R.D- 
OT  IS   PI  TOM  RRTE     UERSUS  TINE 


6 


15? 


PLOT   51 


010  015  020  025  030 


iCP.LE--;5.00E+00  UNITS  INCH.  153 

iCRLF.--=l-OOE+00  UNITS  INCH- 

iOGRRh  3,  TURN  RT  20  KNOTS  WITH  R*D»  6 
-OT  "IS  ROLL  RNGLE       UERSUS  TINE 


1X3 


PLOT   5? 


3CnLF.-5.00F_+00  UNITS   INCH. 
SCnLE>5.00E>00   UNITS    INCH, 


3QGRRM 
-0T  15 


Q 


TURN  m  20  KNOTS  WITH  R.[ 
ROLL  RRTE     UERSUS  TIME 


6 


154 


PLOT   53 


0J.0 


CIS 


02Q 


C2E> 


030 


5CRL5>ro.OGr.>00   UNITS    INCH., 
CRLfc>2.00E-01   UNITS    INCH, 

IQGRRM   3,    TURN   RT   20   KNOTS   NITH   R»D 

UERSUS   TIME 


7 


OT  is       rch     NOLI 


155 


PLOT  54 


3CRLE=5.00E+Q0  UNITS  INCH. 
3CRLE=1.00E+00  UNITS  INCH. 

}0GRRM  3,  TURN  AT  20  KNOTS: WITH  R*D 
-QT  IS   PITCH  RHTE     UERSUS  TIME 


156 


plot  55 


010 


020 


:25 


030 


. 


CRLE---5.0QE+0G  UNITS  INCH. 
CRLF>1.QQE+Q0  UNITS  INCH. 

10GRRM   3*    TURN   RT   20   KNOTS   WITH   R*D 
-QT    IS   ROLL   ANGLE  UERSUS   TIME 

157 


7 


PLOT   56 


030 


3CRLE-5.00E+00  UNITS    INCH. 
3CRLE-5.00E+00  UNITS    INCH. 

10GRRM  3,  TURN  RT  20 
OT  IS    ROLL  RRTE 


KNOTS  WITH  R.D. 
UERSUS  TIME 


7 


158 


PLOT   57 


010 


ij-O 


ozo 


025 


05Q 


3CHLE>5.00E+00   UNITS   INCH. 
3CRLE->2.00E-01  UNITS   INCH. 

}QGRRM   3.    TURN   RT   20   KNOTS   WITH   R.D 
-OT    IS   PITCH   ANGLE  UER5US   TIME 


8 


159 


PLOT   58 


iCRLE=5.00E+00  UNITS   INCH. 
iCRLE-l.OOE+OO  UNITS   INCH. 

QGRflM  3,  TURN  AT  20  KNOTS  NITH  R.D. 
OT  IS   PITCH  RATE     UERSUS  TIME 


8 


160 


plot  59 


OiO 


020 


025 


030 


- 


SCALED.  00E+00  UNITS    INCH 
3CRLE=1.00E+00  UNITS    INCH 

IQGRRM   3,    TURN   RT 
-QT    IS   ROLL   RNGLE 


161 


0   KNOTS   WITH   R^D^    8 
UERSUS   TIME 


n>\ 


PLOT  60 


A 


\ 


\7\J    hf    V    W^^^02C C2T 


030 


)CRLE>--5.00E+00  UNITS    INCH. 
)CRLE>5.00E+00  UNITS   INCH. 

iQGRRM   3o    TURN   RT   20   KNOTS   WITH   R»D»    8 
QT    IS  ROLL   RATE  ULRSUS   TIME 
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